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ABSTRACT 


The application of vector notation to generalized flutter equa- 
tions is described, and the advantages of graphic stability solu- 
The expressions for three-dimensional 
The use of 


tions are pointed out. 
aerodynamic and mechanical coefficients are given. 
stability curves in the presentation of data and the interpretation 
of these curves in terms of flutter characteristics are explained. 

For solving the determinantal equation analytically, a con- 
venient method of successive reduction in order is presented. 
In graphic solutions certain circle relationships have been found 
useful; when applied to determinantal elements, these relation- 
ships form the basis of solutions developed in two and three de- 
grees of freedom. 

The analysis of a symmetric flutter mode of a twin-fin tail is 
outlined to illustrate the methods of solution described. Actual 
examples of the evaluation of parameters are given and a circle 
diagram is applied to the solution of the stability equation. The 
circle diagram and a two-degree-of-freedom M-curve are used to 
determine the effects upon flutter characteristics of variations in 
parameters. 

The simplifying assumptions as applied to different problems 
vary widely. Likewise, variations in the detailed graphic opera- 
tions are used, depending upon the nature of the data and the 
information desired. 


INTRODUCTION 


p | ‘HERE HAS previously been developed a flutter sta- 
bility equation in generalized coordinates which is 


applicable to any number of degrees of freedom.! 
Occasionally, it is desirable to solve this equation for 
as many as four or five degrees of freedom, but for most 
routine problems only two or three degrees of freedom 
need be considered simultaneously. As might be ex- 
pected, the labor required in setting up and in solving 
the equation increases rapidly with the order of the 
determinant. 

Because of assumptions made or because of lack of in- 
formation, some of the flutter parameters are often 
somewhat uncertain. Other parameters actually vary 
under different flight conditions. It is important, then, 
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that the methods of solution employed permit evalua 
tion of errors and estimation of the effects of changes in 
parameters. ; 

In the calculation of the stability determinant con- 
siderable labor is saved by expressing the complex ele 
ments in terms of magnitude and angle rather than as 
With the aid of a few 
relationships 


real and imaginary numbers. 
simple geometric 


handled by graphic methods whenever possible. 


these vectors are 
Thus, 
well-known advantages in making rapid approximate 
calculations, in eliminating large errors, in estimating 
accuracy, and in providing the analyst with a feeling 
for general effects are utilized. Similar methods have 
been applied to specific types of problems by Kassner 
and Fingado,? Bergen and Arnold,* and Arnold.‘ Ex- 
cept in the last case, however, these methods do not 
appear suitable for three-dimensional problems, and 
in the last case a change in the value of a parameter 
would seem to require almost a complete repetition of 
the solution. 

The vector methods of analysis presented here have 
been made applicable to as large a variety of problems 
as possible. Based upon a single set of aerodynamic 
coefficient plots, they eliminate virtually all other 
preliminary work, such as development of specialized 
charts. Most important, the methods have been de- 
signed to meet the requirements mentioned above. 
They are applicable to any number of degrees of free 
dom, they allow ready estimation of errors, and they 
permit rapid evaluation of the effects of changes in 
parameters. 


SYMBOLS 


= aerodynamic force coefficient 
airfoil half-chord 

= reference half-chord 
distance from airfoil mid-chord to control-surface 
hinge, as a fraction of b 

distance from airfoil quarter-chord line to control- 
surface hinge 


diameter of circle expressed as a complex vector 
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F = structural force coefficient 
normalized deflection mode 


= structural damping coefficient expressed as a fraction 


of stiffness 

h = value of f; at quarter-chord line 

a = imaginary part of vector V;; 

1 =V+* cm | 

4,j = general subscripts indicating deflection modes 

k = reduced frequency (= bw/v 

ko = reduced frequency referred to 

M;; = mechanical force coefficient 

mi; = inertia term 

n = number of degrees of freedom 

R;, = real part of vector Vj, 

r = distance aft of control surface hinge 

ro = radius of gyration of total airfoil section about its 
center of gravity 

rp = radius of gyration of control-surface section about its 
hinge 

S = reference span 

Vi; = vector in the complex plane 

Vp = product vector 

t = true air speed (= u0/k+1/X 

vo = reference speed (= woo 

X = flutter frequency parameter (= «7/w? 

x = chordwise coordinate measured aft from quarter-chord 

Ss = x-coordinate of total airfoil section shear center or 
elastic axis 

Xo = x-coordinate of total airfoil section center of gravity 

xg = distance from hinge to center of gravity of control- 
surface section 

y = spanwise coordinate 

= spanwise coordinate of vertical surfaces 

a,, 8; = weighting factors 

p = air density 

po = air density at sea level 

6; = angle of vector Vj; 

Q = structural coupling term 

2 = primitive frequency parameter (= w;?/ wo?) 

w = circular flutter frequency 

@; = natural or primitive frequency of the mode f; 

wo = reference frequency 

Ay = element of span of airfoil 

Am = mass of airfoil segment Ay 

Amg = mass of control-surface segment Ay 


VECTORIAL DETERMINANT ELEMENTS 

The general equation for flutter stability is expressed! 
by the following determinant: 
Ay + My 
d { 90 + Mo 


A In + Mi, 
Ax, + Me, 


Au + Mu 
An + Ma 


An + Mani Aj». + Mi Bes + Man| 


The three-dimensional aerodynamic coefficients in 
this equation are given by the general expression 
. / ” / A 
Dla;’a;"Aag + ai'h;"Aan + a:'B;"Aag + 


h,'c;"Aea + hy'hj"Acn + hi'B;"Aca + 
B,' a;" Ava + B,'h,;"A bh + B,'B;"A bal ( Ay S) (2) 


Ay = 


Since this paper will be concerned with only a few of the 
more common flutter examples, the tab functions have 
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been omitted from Eq. (2). The drag term, being of 


little importance, has also been omitted. 





The two-dimensional A’s in Eq. (2) are Theodorsen’s 
In the 
present development the choice of a value of Theo 


coefficients® with the mechanical terms omitted. 


dorsen’s a in these expressions is completely arbi 
trary and amounts to choosing a reference axis for de 
fining moment coefficients and normal displacements. 
Some of the aerodynamic coefficients, notably A, 
and A,a, assume particularly simple forms if the 





quarter-chord line is used for this purpose; accordingly, 


a is put equal to —'/s, and the following expressions 
are obtained for the real and imaginary parts: 





| tid fe : ae Fe 
vg 1 k?\ 4 T 
Tag = (1/k)(—(2p/x) — (T;/7)] 
Ra = — */s In. = 0 
T; T; 1 Tio, T 2 F 
RR. = (c i 2) -G = 
T a 8 7 kw T w k? 
hi ] Te py Tale ( = i= ry | 
kL a k T 1 2r/ | 
T3 Ll (7; a) TT 1, 
R, = =“ =— 14 
1 t k? \ x? T° 27? 5° 43 
TieT 9 F 
Ww? =k? 
] Tel , T2T 10 l . TT 
i, = f G 
, k| 2x? bd r> ok 27? 
R) —= (7, 7. =— {(Tie w)(1 k)G| 
Ton = (1 k)(T 2 w)F 
Rea = —'/2 — 2(1/R)G + 2F/kR? 
Tea = (1/k)(2F + 2(1/R)G + 1] 
T Ti, 1 Ty F 
ha SS oe SG ete 
4 nx k * nw k? 
ip oo T; 
I = f 2 63 G a 
e kL ¥ nr k T 
Rea —1— 2(1/k)G In = (1/k)2F 


The parameters in these equations have been defined 





and tabulated by Theodorsen and Garrick.’ Numerical j d 
values of the coefficients are given in Tables 1 to8. It e 
is again emphasized that coefficients based on any t] 
other reference axis besides the quarter-chord line are p 
obtained from these in accordance with Eq. (2). ec 

In many airplane flutter problems, particularly hi 
when the taper of the surface over the span to be con- fa 
sidered is slight and an effective rectangular planform hi 
can be assumed, the two-dimensional coefficients ap- b) 
pear singly or in simple combinations in the elements ie 
of the determinant. For these problems, the aero- ea 
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TABLE 1 
Two-Dimensional Aerodynamic Coefficients Independent of c 


Raa Rah Tah Rea 


—0.50 0 —U.9 


—0.50 O —(). 186 


—0.50 0 0.382 


GW GW Ww 
Qu 
or or ¢ 


$4 
64 


50 0 1.523 


—0.50 0 3.175 
—0.50 ) $+. 886 
8.138 

11.714 
15.473 
25.319 
37 . 766 
61.437 

114.49 

169.35 

272.41 

16.667 $99 8&5 


36 


3 
25 
16 
1156 


09 


WwWwwwww w 


— oo os on aN SS SS SY Sy a Se 


or or or 


0.24 0576 


0.20 O4 


0.16 0256 
0.12 O144 
0.10 010 
0.08 0064 
0.06 0036 
0.05 0025 
0.04 0016 
0.025 000625 
0.01 0001 
0 


20.0 5 731.92 


25.0 5 1,163.6 


1s] 


or or or or 


40.0 3,060.9 
100.0 50 19,657 — 667.52 


TABLE 2 


Two-Dimensional Aerodynamic Coefficients, ¢ = 0 


Rag Tag Rba Iba Rog Ing Roh 
—0.169 —0.169 —0.082 —0.212 
—0.089 . 356 —0.147 .211 —0.051 160 —0.208 
+0.052 . 594 —0.108 .349 +0.005 264 —0.202 060 
0.329 890 —0.030 .515 +0.115 . 889 —0.192 095 
716 . 187 +0.082 3 274 . 506 —0.181 132 
105 424 0.199 .793 .437 594 —0.171 163 
821 781 0.421 . 963 744 .716 —0.156 213 
585 O77 815 —0.142 260 
3.368 425 897 —0.131 303 1: 
358 326 049 —0.106 398 20 
7.789 497 30.630 
12.265 
21.936 


31.662 


095 0.666 102 
374 0.923 218 
968 595 440 
561 2.445 633 
451 062 .876 
935 687 2.177 
122 . 434 2.344 
49. 567 903 8.474 2.497 
88. 250 870 34.024 2.555 
127.155 244 863 2.474 63.603 259 .144 


198.775 17.805 9.351 2.163 100. 57 0.797 184 


‘ia 


528 


452 .176 —0.083 
570 Reg —0.052 651 19 _ 886 
259 .491 007 .918 93.137 


rwwo~_ = 


Wwe O11 QO OG 


or WW wee me 


062 561 023 137 137.92 


025 . 585 062 469 222.11 


673 462 112 031 108.04 
{84 597.83 
165 951.56 
216 2,503.1 


422 16,083.5 


Se 


_ 
— 


509. 127 .488 96 399 262.26 —1.268 264 
3,182.9 71.220 342.7 818 1,669.4 845 0.446 


x x 


Go Gr GW bo bo 


oO x —« « —« 


dynamic coefficients are permanently plotted on a scale coefficients thus revised are designated Aq’, Ayn’, and 
convenient for measuring vectors. Typical curves of Ap’. 

this type are shown in Fig. 1. Values of the inde- When taper is to be taken into account, the parame- 
pendent variable 1/k are shown along the curves. All ter k varies over the span in proportion to the chord. 
coefficients relating to control surfaces vary with c, the The summation in Eq. 2 is then carried out separately 
hinge location parameter, and so are represented by for the real and imaginary components, and the result- 
families of curves, as in Fig. lb. The control-surface ing three-dimensional coefficients are plotted and used 
hinge-moment coefficients (those having the subscript in the same manner as the two-dimensional functions. 
b) are generally multiplied by 2/(1 — c)* for conven- The mechanical coefficients of the stability deter- 
ience in plotting, since this factor brings the curves of minant are defined! as 

each family to the same order of magnitude. The My = —my({1 — 25X11 + 2g;;)] (3) 
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TABLE 3 
Two-Dimensional Aerodynamic Coefficients, c = 0.1 
1/k Rag Tag Roba lha Rog Iba Ron Ton Rep lis 
0 —0.135 0 —(0.135 0 —0.054 0 —0.165 0 —0.165 0 
0.5 —0.048 0.323 —0.118 0.165 —0.029 0.118 —0.162 0.026 0.064 0.446 
0.833 +0.107 0. 538 —0.089 0.272 +0.016 0.195 —0.158 0.046 0.484 0.699 
1.25 +0.410 0.807 —0.030 0.402 0.106 0.286 —(0.150 0.072 1.337 0.936 
1.667 0.833 1.076 +0.055 0.525 0.234 0.373 —0.142 0.100 2.583 1.068 
2.0 1.259 1.291 0.144 0.619 0.366 0.438 —0.134 0.124 3.883 1.094 
3.5 2.043 1.614 0.312 0.753 0.613 0.528 —0.123 0.162 6.369 0.998 
2.941 2.879 1.898 0.497 0.862 0.882 0.593 —0.113 0.197 9.117 0.782 
3.333 3.736 2.151 0.691 0.954 1.161 0.660 —0.104 0.229 12.016 0.486 
4.167 5.914 2.689 1.220 1.130 1.885 0.772 —0.085 0.301 19.638 —0.414 
5.0 8.575 3.227 1.843 1.284 2.788 0. 866 —0.068 0.376 29 . 307 — 1.667 
6.25 13.474 4.034 3.066 1.480 +. 482 0.978 —0.044 0.493 47.749 —4.120 
8.333 24.059 5.379 5.808 1.727 8.227 1.100 —0.010 0.695 89.193 —9.438 
10.0 34.704 6.454 8.643 1.870 12.057 1.153 0.013 0. 860 132.12 — 14.558 
12.5 54.301 8.068 13.969 2.010 19.196 1.173 0.042 L712 212.83 —23.385 
16.667 96.641 10.757 25.992 2.094 34.829 1.088 0.080 1.537 391.11 — 40.384 
20 139.221 12.909 37.716 2.066 50.672 0.944 0.104 1.879 573.11 — 55.459 
25 217.609 16.136 60.026 1.878 80.042 0.614 0.134 2.395 912.36 —80.272 
40 557 . 289 25.818 158.075 0.689 208.41 —0.870 0.195 3.946 2,400.4 — 164.27 
100 3,483.6 64.544 1,015.76 —11.756 1,324.6 —12.839 0.333 10.154 15,425 —617.03 
x co x — cx x — =x x x —_ 
° (a) 
v i) ] 
A oe 
ax! /A 
- cha Ace mr 
“4 
8 
Sj ao 
IK 20 ae if 
/ 7 9 Agh?-050 RQ + 
/ 4 a 
9° al xt 
Mal £ 
] 21 / 10 P 
an 
| i 
{ef el Ee ee Le er oe Soe 0 | 
-1,0 O 1.0 20 SD 6 5.0 60 - 
10) rr Mog" 
Fic. 2. Vectorial determination of a determinant element 
When vector notation is employed, any determinant 
element A;; + W/;; is expressed by a single vector and 
can be found in a single operation through use of a scale 
and protractor. As an example, Fig. 2 shows how the 
. ° 7. 
element Aq + M,q is measured; the operation is the 
10 same for the other coefficients. 
- “ta 
THE STABILITY CURVE 
Fic. 1. Typical two-dimensional aerodynamic coefficients 
Ihe vectors representing the determinant elements 
vary with both ky and X. For a given value of kp the 
where tatethtte of : ; a . oe 
stability determinant becomes a complex equation de- 
1 f [x02 + 1? X04 a eN termining two quantities, one of which is the unknown 
n —OjQ; yh 1;a2;) - , oP os . 
me T pby*S l i “er _ X. For the other unknown it is convenient to choose 
; [ets + rg? one of the structural damping coefficients Gi. Both X 
hh; | Am + bo? (a8; + Byay) + and g; can then be considered as functions of fp. 
0 ‘ e f i = ‘ 
. +22 ) The curve showing g; as a function of 1/koV_X or ol 
+8 < e ° 66 ons cle tan 
j (Bh; + hiB;) + ‘c 2.3, | Ams f (4) the true air speed is called the ‘‘stability curve’; it ind1- 
0 o~ . le 
cates the damping necessary to maintain steady-state 


Thus, when there is no damping, the mechanical terms 


are entirely real. 




















oscillations at any speed. If a horizontal line is drawn 


to indicate the true damping coefficient of the structure 
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TABLE 4 


Two-Dimensional Aerodynamic Coefficients, c = 0.2 














FLUTTER ANALYSIS 


443 


Icp 








line and the stability curve give the critical flutter 
speeds. Since the structural damping generally in- 
creases with amplitude, the height of the curve above 
the true damping line is a measure of the violence of 
flutter; where this height is negative, its magnitude 
can be interpreted as a ‘“‘margin of stability.”’ 

The use and interpretation of the stability curve de- 
pend upon a proper choice of the variable g,;; the one 
selected should be that most effective in preventing the 
mode of flutter under investigation. When there are 
several possible modes, it may be advisable to choose a 
different variable g for each mode. For instance, there 


1/k Rag Tag Roba Iba Rog Tog Roh Tbh Ree 
0 —0.104 0 —0.104 0 —0.034 0 —0.125 0 —0.125 0 
0.5 —0.011 0. 287 —0.092 0.124 —0.015 0.084 —0.123 0.019 0.091 0.376 
0.833 +0. 156 0.478 —0.071 0.205 +0.021 0.138 —0.119 0.033 0.486 0. 584 
1.25 0.481 0.716 —0.027 0.304 0.092 0.203 —0.114 0.053 1.290 0.765 
1.667 0.935 0.955 +0.035 0.397 0.193 0.265 —0.107 0.073 2.469 0.839 
2:0 1.393 1.146 0.100 0.469 0.296 0.311 —0.102 0.091 3.700 0.830 
2.5 2.235 1.433 0.223 0.570 0.490 0.375 —0.094 0.119 6.057 0.677 
2.941 3.133 1.686 0.359 0.653 0.701 0.426 —0.086 0.144 8.665 0.417 
3.333 4.054 1.910 0.502 0.723 0.919 0.469 —0.080 0.168 11.418 0.091 
4.167 6.393 2.388 0.875 0.858 1.484 0.548 —0.066 0.221 18.660 —(). 883 
5.0 9.252 2.866 1.347 0.978 2.188 0.615 —0.053 0.276 27 . 856 —2.190 
6.25 14.516 3.582 2.245 1.169 3.506 0.694 —0.036 0.362 45.402 —4.702 
8.333 25.886 4.776 4.257 1.326 6.413 0.781 —0.011 0.510 84.855 — 10.065 
10.0 37.322 5.731 6.338 1.443 9.382 0.819 +0.006 0.631 125.72 — 15.187 
12.5 58.374 7.164 10.247 1. 563 14.909 0.835 0.028 0.816 202.61 — 23 . 963 
16.667 103. 857 9.552 18.874 1.654 26.992 0.779 0.055 1.128 372.43 —40.777 
20 149.600 11.463 27.675 1.656 39.231 0.681 0.073 1.379 545.81 — 55.639 
25 233 . 809 14.328 44.049 1.555 61.924 0.453 0.095 1.758 868.43 —80.027 
40 598.716 22.925 116.01 0.787 161.00 —0.574 0.140 2.896 2,286.7 — 162.30 
100 3,742.5 57.313 745.50 —7.925 1,021.6 —8.903 0.241 7.452 14,697 — 602.67 
co oD co @ — co @ — © oD co oo — = 
TABLE 5 
Two-Dimensional Aerodynamic Coefficients, c = 0.3 
1/k Rag Tap Roa Iba Rog Tog Roh Toh Rep Icp 
0 —0.077 0 —0.077 0 —0.020 0 —0.091 0 —0.091 0 
0.5 0.021 0.247 —0.069 0.091 —0.005 0.057 —0.089 0.014 +0.110 0.309 
0.833 0.197 0.412 —0.054 0.150 +0.022 0.094 —0.087 0.024 0.479 0.473 
1.25 0.539 0.618 —0.023 0.222 +0.076 0.138 —0.083 0.037 1.232 0.602 
1.667 1.019 0.823 +0.021 0.291 0.154 0.179 —0.078 0.052 2.338 0.633 
2.0 1.502 0.988 +0.067 0.343 0.233 0.210 —0.074 0.064 3.494 0.582 
2.5 2.390 1.235 0.154 0.417 0.380 0.253 —0.069 0.084 5.710 0.381 
2.941 3.337 1.453 0.250 0.479 0.540 0.288 —0.063 0.102 8.164 0.082 
3.333 4.309 1.647 0.351 0.531 0.706 0.317 —0.059 0.119 10.757 —(0.274 
4.167 6.776 2.058 0.615 0.631 1.134 0.371 —0.050 0.156 17. 584 —1.297 
5.0 9.791 2.470 0.948 0.720 1.666 0.416 —0.040 0.195 26.256 —2.640 
6.25 15.243 3.088 1.582 0.884 2.661 0.469 —0.028 0.255 42.815 —5.178 
8.333 27.335 4.117 3.003 0.986 4.851 0.529 —0.010 0.360 80.065 —10.527 
10.0 39.400 4.940 4.473 1.078 7.083 0.555 0.002 0.446 118.67 —15.597 
12.5 61.601 6.175 7.234 1.178 11.234 0.567 0.017 0.576 191.29 — 24.238 
16.667 109.57 8.234 13.327 1.267 20.299 0.531 0.037 0.797 351.74 —40.711 
20 157.82 9.880 19. 544 1.288 29.481 0.467 0.049 0.974 515.56 — 55.222 
25 246. 64 12.351 31.109 1.246 46.474 0.317 0.065 1.241 820.39 —78.973 
40 631.51 19.761 81.9386 0.792 120.61 —0.364 0.096 2.046 2,160.7 — 158.85 
100 3,947.3 49.402 526. 55 —5.007 764.21 —5.802 0.168 5.264 13,887 — 583.75 
oO ox oc -_ x — © x x x — 
for low-amplitude oscillation, the intersections of this are two common flutter modes in wings, one having a 


much lower frequency than the other. In the low-fre- 
quency mode the damping term in bending has the 
largest stabilizing effect, while in the high-frequency 
mode torsional damping is generally most effective. 
Hence, in this case it is convenient to solve for ge 
(bending) for the larger values of X, and g; (torsion) 
for the smaller values. The two curves may have an 
overlapping flutter range. 

If Eq. (1) were solved for any M,; and only the real 
parts of each side retained, g; would be eliminated and 
an equation in X of degree 2NV — 1 would result, where 
N is the number of diagonal terms containing X. 
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1/k Rag lop 
0 —0.055 0 
0.5 0.048 0.208 
0.833 0.229 0.346 
1.25 0.584 0.520 
1.667 1. O80 0.693 
2.0 1.579 0.831 
2.5 2.498 1.039 
2.941 3.476 1.233 
3.300 4.484 1.386 
4.167 7.036 1.732 
5.0 10.156 2.078 
6.25 15.900 2.598 
8.333 28.309 3.464 
10.0 40.789 4.157 
12.5 63.764 5.196 
16.667 113.41 6.928 
20) 163.33 8.314 
os 255.23 10.392 
10 193.45 16.628 
100 4,084.3 41.569 
1/k Rag las 
0 —(0).036 0 
0.5 0.068 0.167 
0.833 0.251 0.278 
1.35 0.610 0.417 
1.667 1.1138 0.556 
2.0 1.618 0. 667 
2.5 2.549 0.833 
2.941 3.541 0.980 
3.333 $. 559 ee 
+. 167 7.1438 1.389 
9.0 10.302 1.667 
6.25 16.116 2.083 
8.333 28.679 2.778 
10.0 41.314 3.333 
12.5 64.573 +. 167 
16.667 114.82 5.555 
20 165.36 6. 666 
25 258.40 8.333 
40 661. 56 13.333 
100 4,135.0 33.332 
co Cc ao 


Hence, for a given value of k) there is always an odd 
number of real solutions for X, but never more than 


2N —1. 
hence, of no physical significance. 


Some of these solutions may be negative and, 
Others may have 


large negative g; values; these solutions are of no prac- 
tical importance because they correspond to conditions 


of motion under large mechanical driving forces. 


In cases of three or more degrees of freedom it is 
generally unnecessary, therefore, to examine all the 
In most cases it is possible from 


roots of Eq. (1). 


physical considerations and past experience to deter- 
mine roughly the frequency and speed ranges that are 
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ABLE 6 
Two-Dimensional Aerodynamic Coefficients, c = 0.4 
Rie Dies Robe Ibs Roh Tp) Res leg 
—0.055 0 —() 01107 0 —0.062 0 —(). 062 0 
—() 049 0,062 0.00015 0.0386 —(0.061 0.009 0.123 Q. 2+¢ 
—0.039 0.103 0.02030 0.059 —0.060 0.016 0.464 0.3¢ 
—(0.018 0.152 0.06023 0.088 —0.057 0.025 1.160 0.449 
+0.012 0.199 0.11703 0.114 —0.054 0.035 2.186 0. 43¢ 
0.042 0.235 0.175 0.134 —0.052 0.043 3.261 0.354 
0.101 0.286 0.283 0.161 —0.048 0.056 5.323 0.113 
0.165 0.329 0.400 0.183 —(0.044 0.068 7.609 —(). 21¢ 
0.233 0.364 0,521 0.202 —0.041 0.080 10.025 —0.59 
0.410 0.434 0.833 0.236 —0.034 0.105 16.392 — 1.647 
0.633 0.496 1.219 0.264 —0.028 0.131 24.485 —3.003 
1.059 0.576 1.941 0.298 —().020 0.171 39.946 —5.5: 
3.018 0.682 3.526 0.337 —0.00840 0.242 74.746 —10.79 
2.999 0.748 5.139 0.354 —(0.00042 0.299 110.82 — 15.748 
1. 852 0.821 8.134 0.362 0.00972 0.387 178.70 —24.152 
8.942 0.891 14.666 0.340 0.02308 0.535 328.69 Bh) 1¢ 
13.114 0.914 21.277 0.301 0.03148 0.654 $81.85 — 54.108 
20.875 0.897 33.503 0.208 0.042 0.833 766.86 — 76.978 
54.989 0.630 86.805 —0.216 0.063 1.373 2,020.0 —153.¢ 
353.389 —3.105 549.06 —3.677 0.111 5.533 12,985 — 558. 62 
TABLE 7 
['wo-Dimensional Aerodynamic Coefficients, c = 0.5 
Rie Iba Rog Tog Roh Loh Res leg 
—0. 036 0 —(). 00539 0 —0.04008 0 —(0. O04 ) 
—().032 0.040 0.00246 0.021 —(), 03943 0.00575 0.128 ).18 
—0.026 0.066 0.01655 0.035 —(0. 03844 0.00994 0.438 0.2 
—0.013 0.098 0.04440 0.051 —0.03681 0.01557 1.075 ) ) 
+(). OO5 0.128 0. O84 0. 067 —(). 038491 0.02174 2.013 259 
0.025 0.151 0.124 0.078 —(). 03330 0.02610 2 998 ) 
0.061 0. 184 0.199 0.094 —(0, 03080 0.08515 SSS —(0).12 
0.102 0.212 0.280 0.107 —(0.02859 0.04279 6. 987 —().469 
0. 144 0.235 0. 364 0.118 —(). 02664 0. 04986 9 POT —(). 858 
0.255 0.280 0.578 0.1388 —(0). 02262 0. 066 15.059 —1.919 
0.395 0.320 0.843 0.154 —0.01887 0.082 22.505 —3.2¢ 
0.661 0.373 1.340 0.174 —().91370 0.107 36.737 —5.732 
1.258 0.441 2.424 0.197 —0.00631 0.151 68.783 —10.827 
1.875 0.489 3.526 0. 207 —0.00131 0.187 102.01 —15.591 
3.034 0.540 5.571 0.212 +-0, 00503 0.242 164.55 — 23.635 
5.593 0.592 10.023 0.200 0.01339 0.335 302.76 —358. 855 
8.204 0.613 14.524 0.178 0.01864 0.409 443.90 — 52.153 
13.059 0.613 22 . 842 0.115 0.02517 0.521 706. 56 —73.838 
34.405 0.476 59.078 —0.117 0.03840 0.859 1861.5 — 146.35 
221.114 —1.743 373.01 —2.105 0.06837 2.210 11,968 .0 — 526.91 
co — © ce — c co co fos] — 


most susceptible to flutter. The corresponding ranges 
of kj and X can then be examined by the methods to be 
described. Thus, all those solutions that are unim- 


portant from a flutter standpoint are easily eliminated. 


SOLUTION OF THE DETERMINANT 


A convenient method of solving the determinant is 
that of successive reduction in order, modified from 
Frazer, Duncan, and Collar.? The scheme is as follows: 

(1) Choose the unknown element. 

(2) Substitute a trial value of X in all elements 
except the unknown. 
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TABLE 8 
['wo-Dimensional Aerodynamic Coefficients, ¢ = 0.6 
L/h Ras las Roa Toa Rog log R Res 1 cp 
0 —(0).021 ) 0.021 0 —(0,.002 0 0.02322 0 0.025 0 
0.5 0.081 0.126 0.019 0.0238 0.003 0.011 —(). 02285 0. OOB2E 0.126 0.132 
0.833 0.262 0.209 —0.016 0.038 0.012 0.018 —(0). 02229 0. 00562 0.40: 1S3 
1.25 0.615 0.314 —0.008 0.057 0.030 0.027 —(0.02137 0.00881 0.972 0.184 
1.667 1 0.418 1) 002 0.074 0.055 0.034 —(). 02030 0.01227 1.813 0.104 
1 0.502 0.013 0. O88 0.081 0.040 -0.019389 0.01522 ? 696 0.022 
o 5 2. 52! 0.628 0.034 0.107 0.129 0.049 —().01797 0.01988 $ 395 0.313 
2.941 3.503 0.738 0.056 0.123 0.181 0.055 —().01672 0.02421 6. 282 0.668 
3.333 1. 506 0.837 0.080 0.137 0.234 0.061 —(0.01562 0. O2820 &. 280 1.056 
4.167 7.052 1.046 0.143 0.163 0.371 0.071 —(0.01335 0.08704 13.550 2 096 
5.0 10. 164 1.255 0.2a2 0. 187 0.540 0.080 —(0.01123 0.04630 20. 260 — 3.391 
6.25 15.895 1. 569 0.373 0.218 0.853 0.090 —(0). 00831 0.061 33.091 —5.750 
8.333 28.273 2.092 0.710 0). 261 1. 536 0.102 —0.00413 0.086 61.997 — 10.570 
10.0 40.723 2.510 1.059 0.288 9 234 0.107 —0._ 00130 0.106 91.979 —15.048 
12.5 63.642 3.138 1.714 0.319 3.823 0.110 0.00228 0.137 148.42 22.581 
16.667 113.16 4.184 3.161 0.354 6. 32 0.104 0.00701 0.189 273.17 — 36.760 
20 162.96 5.021 4.637 0.369 9.15% 0.093 0.00998 0.231 4100. 58 49.146 
25 254.63 6.276 7.382 0.375 14.376 0.066 0.01367 0.295 637.70 69.274 
10 651.88 10.042 19.449 0.314 37.110 —(). 057 0.02115 0.486 1,680. 4 — 136.37 
100 4074.4 25.104 125.00 —(). 872 233.87 1.069 ) O3809 1.250 0.803 186.55 
(3) Read off the magnitudes and angles of these ele- or 
ments from aerodynamic curves. ul | V,/ an 
“ . ‘o~ > ss = 4 H P (4 ) 
(4) Star an element in a different row and column 
from the unknown; the subsequent work is simplified When the vectors — I’p and Vp’ are laid off from the 
if the starred element does not contain X and if itis a point on the A; curve, the closing vector is 1J;. From 
comparatively large element. the relationships 
(5) Determine the new elements of the determinant , 
. . Mi (real) My(4 OX 
of next lower order as illustrated below. @) 
(6) Repeat steps (4) and (5) until the determinant is Mj, (imaginary ) my QyeX 


reduced to second order. 
The following third-order determinant is given as an 


example of this method. Here the element Ay, + Wu, 


or Vy, is taken as the unknown, and J is starred. 
Each element A;; + .W;; is written as a single vector 
y 

Vu Vie V3 

Vo V22 Vo3*| = O 

] 31 ] "39 J "33 


Reduction to second order is as follows: 


. VisVa yg, V13 Ve 
Vn — — Vie = 
il Ves 12 Vos 
= — (5) 
Vas ns Vas Va Pea V'sa V'2e 
V33 V93 
or 
Vn — V; V2" 
= & (6) 
V3’ V30" 


With Vy replaced by An + Mu, the second-order 


determinant may be readily expanded to the form 


Ay + My —_ Vp —= Vio" Vay’ V2" 


= Vp’ 


g, and X can then be found. 

If the calculated X differs from the assumed value, 
the above procedure is repeated as a trial-and-error 
problem until agreement is reached; the damping 
coefficient corresponding to the final value is deter- 


mined. If the unknown is chosen as previously dis- 
cussed, convergence of the solution is _ usually 
rapid. 


Once X is determined, the true air speed for the given 
set of conditions may be calculated from the relation- 
ship 


v = dow /k\VX (9) 


From Eqs. (1) to (4) it may be seen that the only 
effect upon the determinant of a change in altitude 
occurs in m;;._ Altitude may thus be varied by dividing 
all mechanical coefficients by p/ po. 

The application of this scheme is not limited to three 
degrees of freedom. For higher order determinants, 
reduction follows the same procedure. 


THE CIRCLE METHOD 
The reciprocal of any vector V = R + iJ is 


1/V = (R — 11)/(R? + I?) 
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This may be expressed in common coordinates as x + 
ty. Then 
x= R/(R?+ Lf), vy = —I/(R? + I) 
Elimination of R yields 
x?+ y? + (y/D) = 0 (10) 


If I is positive, this is the equation of a circle below 
and tangent to the x-axis at the origin and having a 
diameter 1/J. As R varies the vector 1/V traverses 
this circle; its exact position is determined by the fact 
that its angle is the negative of the angle of V. 

The above geometric relationship has a useful appli- 
cation in flutter work. For example, the quantity 


—1/V2 = —1/[Ax ad Moo(1 — 2X )] 


is a vector that, as X varies, traverses the circle whose 
diameter extends from the origin to the point +1/J2 
on the imaginary axis. Henceforth, this diameter will 
be treated as a vector that defines the circle. In addi- 
tion or multiplication the diameter may then be treated 
as any other vector. 

In the simple problem of binary control-surface 
flutter expressed by the equation 


Vo Voe 


where V2 and V2; do not contain X, if the solution is 
written in the form 


— My = An — (VeVa1)/ V2 (11) 


the right-hand side can be represented by a circle 
having a diameter vector drawn from the point A). 
The diameter will have a magnitude | Vj_V2;|/Jo2 and a 
direction 6j. + 6 + 90 deg. The exact point on the 
circle corresponding to a particular value of V2 is de- 
fined by a line drawn from A, at an angle of 90 deg. — 
422 to the diameter vector, as in Fig. 3. The point so 
defined is thus the solution to Eq. (11) and gives the 
components of JV). 

Information of a general nature can be drawn from 
this circle alone. For example, if the circle does not 
pass below the real axis, flutter cannot occur for any 
primitive-frequency ratio at the particular value of 
1/ko chosen since no positive g; Solutions exist. 





I ihe 




















> 


Fic. 3. Stability solution in two degrees of freedom by means of 
the circle construction. 


If either Qe or m2 is uncertain, as for example in a 
cable-connected control system where cable tension 
and pilot restraint vary during flight, a conservative 
solution may be found if only the lower extremity of the 
circle is used. Eq. (3) applied to /,, then determines 
both X and g:. 

When V2 contains a damping term ge, variations in 
X cause V2 to traverse an oblique instead of a hori 
zontal straight line. The slope of the line is go. In this 
case the diameter vector is rotated through — tan! go, 
and its magnitude is determined by Je instead of Jo», 
where Jo: is the perpendicular distance from the origin 
to the above straight line. 

The circle construction is useful in exploring the 
effects of changes in parameters, particularly when 
both g,; and ge can be put equal to zero temporarily. 
Then both the above circle, which is called the (11) 
circle, and a similar circle (22), obtained by reversing 
the roles of Vj, and V2 in the above process, can be 
constructed as shown in Fig. 4. The real axis inter- 
sections are values of —M,, and — Mz which satisfy 
the determinantal equation. It will be found that the 
left-hand intersection of the (11) circle combines with 
right-hand intersection of the (22) circle to form one 
solution, the other two intersections constituting a 
second solution. These two solutions are plotted on 
My, vs. Mz. coordinates. Curves in this plane are ob 
tained by allowing other parameters, usually kp or Jp, 
to vary. The curve shown in Fig.'5 is for constant ko 
and variable Vo. 





M,2*Ma)= -0.5 
Vk =4.| 




















Fic. 4. Two-circle solution used in construction of /-curves 





A, 74. 
Lp, +0.76 
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Fic. 5. M-curve as used to study effects of variations im 
. mechanical parameters. 
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Since both /;; and Wy are linear in X, combinations 
that are consistent with Eq. (3) lie along a straight line 
passing through the point (#7, mx). Intersections of 
this line, called the pivot line, with the /-curve there- 
fore represent particular solutions to the stability equa- 
Variations in the primitive-frequency ratio cause 
Variations 


tion. 
the line to pivot about the point (71, 22). 
in My OF Mz shift the pivot point. Variations in Wj. 
or ky are shown on the .W-curve. 
rameter X is a linear function of the distance of the 
M-curve intersection from the pivot point and is easily 
Hence, for zero structural 


The frequency pa- 


determined for all solutions. 
damping, all possible solutions obtainable by varying 
the other mechanical parameters can be found. The 
construction of /-curves is remarkably easy, for when 
A, and Az. are known only the diameter vectors must 
be calculated. 

The application of the M-curve is not necessarily 
limited to zero structural damping. However, the 
introduction of a variable gz; or go usually makes it 
desirable to restrict one of the other quantities ordinar- 
ily allowed to vary. 

The circle construction has been found useful in 
three- and four-degree-of-freedom problems. The 
method of successive reduction in order described earlier 
is used, but in this case a diagonal element containing 
X is starred. Then, if the nondiagonal elements oc- 
curring in the same row and column do not contain X, 
or if they vary slowly with respect to X so that approxi- 
mate values can be assigned to them, it is clear that the 
determinant of the next lower order consists of the 
minor of the starred element but with each aero- 
dynamic term therein replaced by a circle. 

In problems of more than two degrees of freedom, for 
each flutter mode there are usually two degrees of 
freedom which are more important than the others. If 
a secondary degree of freedom has a small influence on 
the total motion because of slight coupling action, the 
circles obtained by starring the corresponding diagonal 
element are small and can be treated as perturbations 
on the vectors that comprise the minor of the starred 
element. If the influence is small because of a rela- 
tively high primitive frequency, the perturbations are 
given by small arcs of the circles, the angle of arc being 
determined as in Fig. 3. If the starred element is 
highly uncertain because of unknown stiffness or in- 
ertia parameters, points that, by inspection, appear to 
represent a condition of least stability may be fixed 
on the circles; these points are then used to obtain a 
conservative solution, as in the two-degree-of-freedom 
case. 


ILLUSTRATIVE EXAMPLE 


Analysts who apply the foregoing principles will 
usually develop personal preferences for particular 
methods, depending upon their immediate aims and 
their previous experience. For this reason the pre- 
ceding discussion has been intended merely as an in- 
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troduction to many possibilities inherent in graphic 
methods of stability prediction. 

As an illustration of typical steps in the solution of 
flutter problems, the analysis of a symmetric mode of a 
twin-fin tail supported rigidly at its plane of symmetry 
will be outlined. 

The coordinates used to describe the motion of the 
tail are as follows: (1) twist of the stabilizer about the 
quarter-chord line, the rate of twist being constant 
and the elevator 
stabilizer bending with constant curvature and a node 


remaining untwisted; (2) pure 


at the plane of symmetry, the outboard vertical sur- 
faces being undeformed; and (3) elevator flapping. As 
defined by the equation! 


fi = hie t+ axe/ho + Bir/do 


the weighting factors are 





k, = 0, a = y/S, Bb, = — y/S 

a2 = Bo = 0, Ae = (y/S)*? over the horizontal surface 
a2 = Be = 0, Ae = 22/S over the vertical surfaces : 
a3; = h; = 0, 8; = 1 over the elevator span 


(12) 


where y is the horizontal surface spanwise station, 2 is 
the vertical surface spanwise station measured from the 
plane of the stabilizer, and S, the reference span, is 
arbitrarily chosen as the distance from the plane of 
symmetry to the plane of the vertical fin. For the 
reference half-chord, do, the stabilizer half-chord at the 
fin-stabilizer intersection is used. 

The following weighting factors are introduced for the 
calculation of aerodynamic coefficients: 


ay’ = (b/bo)*(y/S) Bi’ = — (b/bo)?(y/S) | 
he’ = (b/bo)(y/S)? over the horizontal surface 

= (b/by)(22/S) over the vertical surfaces e9) 
83’ = (b/bo)? over the elevator span 


The broken 
curves represent the 
primed weighting factors! and take into consideration 


These are shown as solid curves in Fig. 6. 
reasonable values of double- 
the spanwise pressure effects resulting from departure 
of the airflow from that of the simple two-dimensional 
theory applied to each element of span; because of this 
aerodynamic interaction; the factors extend over both 
the horizontal and vertical surfaces. 

The three-dimensional aerodynamic coefficients are 
now obtained from Eq. (2). 


Au = Dla:’a:”"Aae + a1’Bi"Aag + Br’a,"Ara + 
81’ 81" A pgl (Ay/S) 

Ay = Llanhe”"Agn + Bi'he” Ay] ( Ay/S) 

Aj3 = =[a1’B3"Aag + 81’ B3”A vg] (Ay/S) 


An = >[he!01"A ca + he’B:"A cp] (Ay/S) saa 


Am = >[he'he”A «| (Ay/S) 
Ass = 2[he'Bs"A eg] (Ay/S) | 
An = 2[6s’01"A re +-Bs’Br"A vg] (Ay/S) | 
Az = >[Bs’h2”A val ( Ay, S) 

A33 = =|83'83"A pg] (Ay/S) J 
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TABLE 9 
Calculation of Ao; for 1/ko = 2.80, c = 0.4 
y/S b/bo 1/k Reg Ics Ay/S hy 3 R I 
0 1.394 2.00 3.26 0.354 0.0472 0 1.943 O 0 
0.094 1.393 2 00 3. 26 0.354 0.0945 0.009 1.940 0.005 0.001 
0.189 1.363 2.05 3. 47 0.330 0.0945 0.049 1.860 0.030 0.00: 
0.283 1.322 y ie | 3.68 0.307 0.0945 0.106 1.750 0.065 0.005 
0.378 1.281 2.18 +. 00 0.268 0.0945 0.183 1.639 0.113 0.008 
0.472 1.242 2.25 $28 0.239 0.0945 0.276 1.542 0.172 0.01 
0.567 1.200 2.32 4.56 0.207 0.0945 0.386 1.444 0.241 0.011 
0.661 1.160 2.40 1.90 0. 167 0.0945 0.508 1.347 0.316 0.011 
0.756 1.119 2.49 5.28 0.122 0.0945 0.640 1.185 0.378 0.009 
0.850 1.080 2.58 5.74 0.058 0.0945 0.780 1.010 0.427 0.004 
0.945 1.0389 2.69 6.24 —0.009 0.0945 0.927 0.825 0.451 —0) 001 
1.040 0.972 2.87 7.20 —0.157 0.0945 1.050 
1.133 0.648 1.30 17.6 —1.83 0.0945 0.802 
z/S 
—0.154 0.411 6.78 17.8 —6.72 0.0236 0.010 —0.039 0 0 
—(0,.094 0.972 2 O7 7.20 —Q0.157 0.0945 0.009 —0.175 —0.001 0 
0 1.162 2.39 1 86 0.172 0.0945 0 +). 390 0 0 
0.094 1.151 2.42 4.99 0.157 0.0945 0.010 0.301 0.001 O 
0.189 1.024 2 .de 6.49 —0.049 0.0945 0.037 0, 223 0.005 0 
0.283 0.802 3.48 11.14 —0.76 0.0945 0.064 0.146 0.010 —0.6001 
0.378 0.440 6.35 41.6 —5.74 0. 0866 0.063 0.049 0.011 —().002 
> 2.224 0.058 
balanced, are assumed to be rigid parts of the vertical 
surfaces. 
The application of Eq. (+) to the calculation of the 
i. inertial parameters is straightforward. For example, 
l xo? + 1%? 
my = — a2 Am + 
a pbo2S by? 
bat. 3 e 2018; + e. 3) Amg (13) 
a bo? by + a 
Moa = am [he? Am| 
rpby*S 
b In the evaluation of these summations, both the hori- 
r zontal and vertical spans are usually divided into seg- 
Osr Ay) ments and the individual quantities computed for each 
OAT B” _4----B3 segment. Sometimes the vertical surface can be treated 
Of Pt Se = as a single unit; in the calculation of my, the moment 
-OAF = . aa — . ‘ x ° 
7 Ree Bi , A of inertia of the vertical surface about a lateral axis 
-O8F eee . ‘a ™ ie 
ie A, can be entered as a single item of the first term. In 
: O2 04 O6 O8 10 2 the calculation of mz. the vertical surface appears 
= i. i. 1 -_ . . - . . 
-0.2 D O02 O04 twice—first, as a single item of weight in the summa- 
We r/c, tion over the stabilizer span and, second, in the sum- 


Fic.6. Variation of aerodynamic weighting factors over the span 
of the horizontal and vertical surfaces of a twin-fin tail. 


The summations extend over both the horizontal and 
vertical surfaces. As an example, the calculation of Ao; 
for a single value of 1/k) is shown in Table 9. Atten- 
tion is called to the fact that A.g varies along the span 
as a result of the taper; also, the aerodynamic coupling 
between stabilizer bending and elevator flapping is in- 
creased by the aerodynamic forces acting on the 
vertical surfaces. The rudders, having a high natural 
frequency for symmetric motion and being fully mass- 


mation over the vertical span. The latter sum, accord- 
ing to the assumed deflection mode, is 4/.S? times 
the moment of inertia of the vertical surface about 
a fore-and-aft axis through the stabilizer-fin intersec 
tion. 

The calculation of the A’s and m’s constitutes a con- 
siderable proportion of the labor involved in any flutter 
problem. However, convenient tabular forms can be 
devised by means of which this work can be assigned 
to computers. 

Primitive frequency calculations are carried out 
without difficulty according to Eqs. (18) and (19) of 
reference 1. In the present example the determination 
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of Qy is particularly simple because the uniformity of 
the the absence of 
make it possible to use the equation 


stabilizer structure and cutouts 


where the elastic axis position x, and the center of 
gravity position x) are determined in the vicinity of the 
fin station. By definition, Q,; 
work done by the system of forces F; Am (real) during a 


Since the system of 


is proportional to the 


virtual displacement equal to f;. 
forces required to deflect the elevator in mode (3) alone 
does no work under a virtual displacement in either 
mode (1) or (2), the remaining structural coupling 
terms 9}3 and Qe3 are zero. 

When all the mechanical parameters are inserted, 
the stability equation for standard atmospheric condi- 
tions at sea level becomes 
Ai; — 0.060 | 


Ax — 0.024 =Q 
A3s—0.117(1—0.50X)| 


Aun—5.51[1— X(1+ig1)] 
Au —2.41(1—0.04X) 
Au— 0.060 


Aw —2.41(1—0.04X) 
Az —5.06(1 —0.20X) 
Az — 0.024 


(16) 


These coefficients have been calculated for an elevator 
that is 33 per cent statically balanced, the balance 
weights being distributed over the outboard third of 
the elevator span. The torsional frequency is taken 
as the reference frequency w. The elevator frequency 
parameter is assumed to be 0.5. 

In the solution of this equation by the circle method, 
it is convenient to star the element V3; and reduce the 
determinant to second order as outlined previously. 
The circles for each of the four reduced elements are 
then drawn (Fig. 7). Since the aerodynamic coefficients 
vary with 1/ky and the mechanical coefficients vary 
with the density ratio, there is a different circle diagram 
for each combination of these parameters. 

For illustrative purposes, the case of 1/po 2.80 
and p/po 0.58 may be considered. The determinant 
is then that of Eq. (16) with all the mechanical pa- 
rameters divided by 0.58. Through use of a scale and 
protractor, the following vector magnitudes and direc- 
tions are determined as shown in Fig. 2: V3; = 0.243, 























150 deg.; Viz = 1.68, 24 deg.; V9 ae 0.067, 160 deg. ; 
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Fic. 7. Circle diagram, showing solution of the vector equation 


in three degrees of freedom. 


OF 
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Vos = 2.18, l deg. 
to be 0.331. 


diameter vectors are 


The imaginary part of V3; is found 
The magnitudes and directions of the 


Vis V3 
Dy, = 7 = 1.23 63 + 0 + 90deg. = 264 deg. 
Vis V3 9 ' - 
Dy = - FT = 0.342 613 + Os. + 90 deg. = 274 deg. 
J 23 V3 . . 
Day 7 = 1.60 Go3 + 03) + YO deg. — 2 }] deg. 
33 
Vos V2 - —_ ~ 
Des = 7 = 0.4 45 Gos = As9 + 90 deg. = 251 deg. 
33 


The diameters are laid off from the respective points 


Ay, Ai, Ao, and Ag, and the four circles are con 
structed. 

The particular arrangement of Fig. 7, with the circles 
drawn on two vertically aligned sets of axes, is useful 
in finding consistent sets of values for the mechanical 
terms Mq;, My, Mo, and Me. All of these coefficients 
are linear functions of X. With —W/, on the upper 
horizontal axis and — My (= —Mz2) on the 
horizontal axis, linearity requires that corresponding 
values of the two coefficients always be on a straight 
Similarly, 


lower 


line passing through a fixed point /P. 
— Me and — M2 must lie on a straight line pivoting 
about a point P,. The two pivots are easily located by 
drawing the lines for any two values of X. For any 
value of —My, the corresponding value of —My 
(= —Mz2) is then given by the intersection with the 
lower axis of the line connecting the points — Wy, and 
P;, and —Mz is given by the intersection with the 


P, 


upper axis of the line through the points and 
— My). 

In the event that My. and 2 are independent of X, 
the circle diagram may be further simplified. The (11) 
circle is drawn on the upper plot; the (22) circle, on 
the lower one; —1/;, and — JA then require but a 
single pivot point. The other two circles may be placed 
wherever convenient. 

To return to the solution itself, it may be seen that 
corresponding points on the four circles are determined 


by the angle of V33, as in Fig. 3. These points are the 


termini of the vectors Vu, Vie, Vor, and V2, the ele 
ments of the reduced determinant, drawn from the 
respective —.W points. If now 

Vii Veo Vie Vo21 (17) 


the determinantal equation is satisfied. But since the 
imaginary part of —J/y is unknown, it is convenient 
to draw Vi from the point —/,, (real) to the vertical 
line through the point on the (11) circle; the length d 
(Fig. 7) then is — Wy (imaginary). 

As a first trial in the solution under consideration, 
it is assumed that — 1, (real) = 0. This value corre- 
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sponds to an X of unity, for which condition 633 is found 
to be 40 deg. The other three M/’s are determined by 
means of the pivot point relationships, the points on 


the circles are located at an angle 90 deg. — 43; = 50 
deg. to the diameters, and the following vectors are 
measured: Vy = 3.98, 182 deg.; Vo, = 1.26, 165 deg.; 


Vee = 6.90, 172 deg. From Eg. (17), Vu is found to be 
0.73, 175 deg. 

The calculated V\,, when drawn from the real axis 
so that its terminus falls on the vertical line R = —0.44 
drawn through the point on the (11) circle, determines a 
new point —J/,, (real) = 0.29. When the trial solu- 
tion is repeated with this assumption, Viz and V2 
show practically no change; V2 and Vy, become 6.95, 
172 deg., and 0.72, 175 deg., respectively; and —My 
(real) retains the above value. The quantities X and g; 
are calculated from Eq. (3) as follows: 


. my + My (real) 9.51 — 0.29 
Vw ce um { ie - = 0.97 


myQh (9.51)(1) 


d —0.21 
= = —0.023 


my + My (real) 9.51 — 0.29 


si. = 


This value of X does not appreciably shift the points on 
the circles as originally located for X = 1.0. 

With the same determinant and the same set of cir- 
cles, it is often possible to find additional solutions 
corresponding to other values of X. The various solu- 
tions can then be approximately examined in order to 
determine the most critical flutter mode. 

The solutions obtained by varying 1/kp and p/po are 
conveniently plotted as shown in Fig. 8. Only one 
flutter mode is represented, since the others are less 
critical in this case. The horizontal scale, 1/ko+/X, 
becomes true air speed when multiplied by the reference 
speed v = dow. Critical speeds for various altitudes 
are determined by the abscissas corresponding to an 
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Fic. 8. Stability curves showing critical speed factors for 
various altitudes and damping factors. 


estimated or experimentally determined value of th« 
structural damping factor g:. 

It will, of course, be understood that high accurac\ 
cannot be obtained by the above process. However 
many parameters are usually so uncertain that an ex 
tremely precise solution is not warranted. If greater 
accuracy is desired, it can be obtained with little extra 
work by direct substitution of the graphically deter 
mined X in the stability equation. 

The flutter characteristics determined from the sta 
bility curve indicate the limits of safe operation for the 
airplane under consideration. If these limits restrict 
performance too greatly, a survey is made to determine 
the design changes that will improve the airplane most 
efficiently. 

The effects of variations in elevator parameters are 
readily observed in the circle diagram, Fig. 7. Changes 
in elevator frequency or moment of inertia merely shift 
the points on the four circles, while changes in balance 
alter the diameter vectors. In this case, however, the 
elevator parameters are found to have little influence 
on the critical flutter speeds. 

For the study of variations in other mechanical 
parameters it is convenient to make Q; infinite and 
thereby reduce the problem to one of two degrees of 
freedom. Suppression of the elevator flapping mode is 
permissible here because of its ineffectiveness in pro- 
moting flutter. Now the /-curve for zero damping 
may be constructed as explained earlier; the sample 
solution of Fig. 4 and the plot of Fig. 5 are actual ap- 
plications from the subject problem. 

The closed curve (Fig. 5) is shown for 1/ky = 4.1, a 
value that is critical for this problem. As an example 
of the use of the curve, pivot lines are shown for 2, = 
0.2 and ® = 0.3, both drawn for p/p. = 0.76. The 
intersections of these lines with the curve determine 
values of M1. for which the stability curves have axis 
intercepts at 1/k) = 4.1. From the shape of these 
latter curves critical speeds can be estimated for par- 
ticular damping factors when the intercepts are 
known. 

From the figure it is clear that for Q; = 0.2 the 
stabilizer unbalance term —W/;. must be reduced from 
3.2 to 1.8 in order to shift the intercept above 1/kp = 
4.1, while for 2, = 0.3 a reduction from 3.2 to 1.2 is 
required. The same result would be accomplished by 
an increase in 2 from 6.7 to 8.1 in the first case and 
from 6.7 to 9.0 in the second, since these changes 
would translate the line QB upward until it intersects 
the curve at — My = 3.2. 

Through use of the M/-curve, then, the effects of any 
changes in structure or mass distribution can be readily 
estimated. For example, added mass along the sta- 
bilizer leading edge would have two important beneficial 
effects—a reduction in — Jy and an increase in 12 
and one minor adverse effect—an increase in 71. 

The fact that the curve is a closed loop shows that 
large values of —W12, as well as small values, are bene- 
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VECTOR METHODS OF 


ficial provided they can be attained without adverse 
effects on the other parameters. 
CONCLUSION 

The type of information desired and the nature of the 
data available vary so widely from one flutter problem 
to another that the methods used must be extremely 
Detailed 
procedures making use of the general ideas presented 


flexible in order to be generally applicable. 


here will occur to the analyst as the need arises in 
particular problems. 

The primary assumptions upon which any problem is 
based are determined by: the nature and accuracy of 
available data; the importance of an accurate solution 
from the standpoint of safety, performance, and manu- 
facturing costs; the status of the problem in relation 
to future designs; the complexity of three-dimensional 
aerodynamic and mechanical effects; and the time and 
man power available. Six degrees of approximation are 
recognized: (1) substitution for the actual system of a 
two-dimensional setup having the properties of a par- 
ticular section of the airfoil—for example, the three- 
quarter span station of a wing; (2) the same, except 
that some or all of the mechanical parameters are calcu- 
lated in (3) calculation of both 
aerodynamic and mechanical parameters in three di- 
mensions with the assumption of no taper or spanwise 
(4) the same, except that 


three dimensions; 


aerodynamic interaction; 
interspan effects such as the tip vortex and interaction 
between horizontal and vertical surfaces are considered 
qualitatively; (5) complete three-dimensional problem 
except that interspan effects are neglected; (6) in- 
clusion of all known three-dimensional effects. 

Most problems involving small, moderate-speed air- 
craft can be treated two-dimensionally by approxima- 
tions (1) and (2). The extensive work of Wylie® * 
makes it unnecessary to do more than determine the 
equivalent two-dimensional mechanical parameters, 
except in cases of more than two degrees of freedom or in 
problems involving tabs. In approximation (3) the two- 
dimensional aerodynamic coefficients appear in differ- 
ent ratios than those in Wylie’s work, and his solutions 
cannot be used. 

For flutter problems that do not readily yield to the 
simplifying assumptions, such as those involving the 
complex modes on extremely large airplanes, one of the 
approximations—(4), (5), or (6)—is used. The fact 
that these solutions generally can be applied qualita- 
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tively to other designs besides the specific model under 
investigation is justification for extended computation. 
It is to be expected that, as such problems are solved 
in detail, short cuts will appear along with a better 
understanding of the relative importance of different 
parameters and that this knowledge can be applied to 
advantage in subsequent analyses. 

In modern high-speed aircraft the effects of com- 
pressibility must be considered. Possio'! has shown 
that in the case of an oscillating wing compressibility 
causes large changes in the aerodynamic coefficients, 
the effects being functions of both speed and frequency. 
In addition, when local speeds exceed the speed of 
sound, the shock fronts and accompanying burbles in- 
troduce nonlinear effects that can best be studied ex- 
perimentally. The existence of these complex aero- 
dynamic effects at extremely high speeds and our 
limited knowledge of three-dimensional airflow make it 
necessary occasionally to supplement theoretic results 
with wind-tunnel tests and flight experience. 
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A Simplified Two-Dimensional Theory 
of Thin Airfoils 


H. J. STEWART* 
California Institute of Technology 


INTRODUCTION 


N THE EXISTING SOLUTIONS of the thin airfoil problem, 
I the determination of the steady state lift and mo- 
ment of an airfoil of small camber and thickness, as 
given by Birnbaum,' von Karman,’ C. B. Millikan,’ 
etc., + > the velocity potential for the flow is determined, 
and the lift distribution and the integrated lift and mo- 
ment are found from this result. It appears that a 
somewhat simpler solution of this problem can be ob- 
tained by working directly with the velocity compon- 
ents rather than with the velocity potential. This same 
method can also be used to advantage in many three- 
dimensional wing problems. An important feature of 
this method is that the inverse problem, the determina- 
tion of an airfoil shape that will produce a given lift 
distribution, can be solved just as easily as the direct 
problem. Both the direct and the inverse problem will 
be investigated by this method in the present paper. 

The notations used in the analysis are as follows: 


x, y = rectangular coordinates in the plane of the airfoil 

r = chord of airfoil 
V -1 

z x + iy 

u, 2 velocity components in the x and y directions, respectively 
u — iv = complex velocity 


complex velocity due to the airfoil 
= complex velocity due to the airfoil camber 


U = velocity of air far from the airfoil 

ae inclination of the flow far from the airfoil 

p static pressure 

Ap differential in static pressure between the lower and upper 
surface of the airfoil 

p fluid density 

( '/o9U? = dynamic pressure 

n airfoil ordinate 

zB Crgc = lift of airfoil per unit span 


Vo = Cyyc? = pitching moment about the center of the airfoil 
per unit span ‘ 
lhe fundamental equations for the two-dimensional 
irrotational motion of an incompressible fluid are 


Ou _ ya O1 - 
Ox oy | 
and 
Ou Ow -. 
Oy Ox = 
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which are the conditions imposed by continuity and ir 
rotationality, respectively. If velocity components are 
found which satisfy these equations for any given 
boundary conditions, the pressure distribution for steady 
state motion may then be calculated from the equation 


b + '/op [u? + v7] = constant (3) 


which is an integral of the equations of motion. If 
Eqs. (1) and (2) are considered as Cauchy-Riemann 
equations, they imply that the complex velocity w = 
“u — iv must be an analytic function of the complex 


variable z = x + 1y, or 
u— iv = w(z) (4) 


If the airfoil ordinates are given by n(x) (Fig. 1), a 
boundary condition for determining w(z) is that the flow 


must not cross the surface of the airfoil; that is 


v = u(O7n/Ox) 5 


on the surface of the airfoil. A second boundary condi 
tion is that far from the airfoil the flow must be a recti 
linear flow of constant velocity. The third boundary 
condition, the Kutta condition, is that the velocity must 
be finite in the neighborhood of the trailing edge. If 
the airfoil is thin and only slightly cambered, the first 
boundary condition may be considered as applying 
along the x axis rather than on the actual airfoil surface. 
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Sketch showing the notation used in the thin 
analysis 


In order to determine w(z), it is convenient to intro- 


duce another complex variable ¢ = re’ where 


The points for which r 1 will be considered as those 


that transform into the z-plane. It is readily seen that 


























TWO-DIMENSIONAL 
points on the circle r = 1 in the ¢-plane correspond to 
points on the x axis in the z-plane where 
x = (c/2) cos 6 (7) 
[he boundary condition of Eq. (5) may thus be con- 
sidered as applying on the unit circle in the ¢-plane. 


LiFT AND MOMENT OF THE AIRFOIL 
rhe total lift and pitching moment of the airfoil are 
most easily calculated by using the Blasius equations. 
If X and Y are the resultant force components in the x 
and y directions, respectively, then 


X — iY = (pi/2) £ w'dz 


where the integration is carried over any contour in the 
fluid enclosing the airfoil. Similarly, the stalling mo- 
ment of the airfoil about the origin is 


M = RI[(p/2) £ wzdz] 


(8) 


(9) 


The velocity field around any airfoil may be expanded 
in a series, valid for large values of £, of the form 


w= U > (a, + ib,) o™ (10) 
n=O 
where the a,,’s and b,’s are all real, ap + ido) = e~ and 
a; = Osince the airfoil is a closed body. From this, the 
resultant lift and moment are 
L = (p/2)U*crb, 
; ’ (11) 
and M = (p/2)U°c*z[ae sin a — be cos a] 
Thus, the lift and pitching moment coefficients are 


C, = 1h, f 
— (12) 
and Cy = 7[a2 sin a — & cos a] 
If the lift distribution is desired, it may be calculated 
fromEq. (3) once the velocity field has been determined. 


THE FLAT PLATE AIRFOIL 


The suggested method will first be used to solve ex- 
actly the problem of the flat plate airfoil. The airfoil 
is placed on the x-axis with its leading and trailing 
2 and x = c/2, respectively, as in Fig. 


edges atx = —c, 


Y 





- 


: ~nil 
a 


NIO 


(3) 





The flat plate airfoil. 


Fic. 2. 
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2. If the airfoil is at an angle of attack a, then w(z) can 


be written as 
w(z) Ue + w’(z) 13) 


where w’ = u’ — iv’ is the velocity due to the airfoil. 


From Eq. (5), v = O and, hence, v = —U sin a = 


constant on the airfoil. w’(s) must also satisfy the 


Kutta condition and vanish as z —~ ©. It is easily 
seen that the function 


(14) 


satisfies the latter two conditions and on the airfoil gives 


op),.( 7) 
wg\e</) 


”» = 1 
Va = Ee 


From this, the flow around a flat plate is given 
by 


ws) = Ue + (21U sin a)/(¢ + 1) 15) 
On the airfoil this gives, of course, no vertical velocity 
and the horizontal velocity is 


u = U[cos a + sin a tan (6/2)| 16) 


If the velocity is expanded in a series of the type used 
in Eq. (10), it is found that for m 2 1, 


a, +16, = 


The lift and moment coefficients for a flat plate airfoil 
are thus 


—27 sin a (—1)’ (17) 


= 27 sin a 


sin a COS a LS) 


- T 
c = 
9 


These are both well-known results. The pressure dis- 


tribution over the plate is given by Eq. (3) as 


p= 


—pU? {sin a cos @ tan (6/2) — 


1/, sin? a [1 — tan?(@/2)}} (19) 
where the pressure has been taken as zero far from the 
airfoil. 


THE AIRFOIL OF ZERO THICKNESS 


Airfoils having zero thickness and a small camber will 
The leading and trailing edges 
+c/2 asin Fig. 3. It will be 


next be investigated. 
are again placed at x = 
assumed that both the-camber n(x) and the angle of 


Y 












Fic. 3. The cambered airfoil of zero thickness. 
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attack a are small quantities of the first order. Squares 
and higher powers of these quantities will be neglected 
as compared to the first-order terms. With this ap- 
proximation and in view of Eqs. (14) and (15), the velo- 
city field may be written as 


w(z) = U(1 — ta) — Uaw) + w’” (20) 


where w” is the velocity due to camber. In the linear- 
ized form the boundary condition of Eq. (5) may be 
written as 

(21) 


” 
v = 


a 


U(0n/0x) 


It can be seen that this linearized boundary condition 
breaks down in the neighborhood of the leading edge 
where x + (c/2) is of the order of ca’. 

Since the airfoil is of zero thickness, v” must be an 
even function of #, and the most general expression for 


w" = Ui > Ag y n (22) 
n=1 


This expression automatically satisfies the Kutta condi- 
tion. This series may not converge on the airfoil, but, 


if it does, the velocities there are 


aa yr |~ vy .° 
u"° =U) )d," sinn 


n=l 


and 


v” = —U > bd," cosn 
n=l 


In order to determine the b,,”’s, a cosine expansion for 


On/Ox must first be obtained. Let 
On/Ox = C)/2 + > c, cos nO (24) 
n=1 
where 
C, = 2/e { @n Ox) cos n6 dé (25) 
0 


As this series contains in general a constant term, ¢o/2, 
the series expansion for w” will not in general converge 
on the airfoil; however it is apparent from the results of 
the analysis of the flat plate airfoil that w” may be 
written as 

(26) 


—[Uie/(—§ + 1)]) — Ui pe n 


n=1 


o 
w= 


where the latter series converges on the airfoil. By ex- 
panding w” in a series of the type of Eq. (10), it is found 


that 


an” + ib,” = i[(—1)"oo — cl 


9; ™ 
= — Ff (2") [cos n@ — (+1)"|d6@ (27) 
a J0 \Ox 


The coefficients for the term — Uawy of Eq. (20) must 
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be added to Eq. (27) to obtain the total lift and moment 


coefficients. These are 


Cr = 27a — 2 [ Oy (1 + cos 6)dé 
J0 \ox 
Cy = x ox +! if (2) (cos 206 — 1)dé 
2 0 \Ox 


These results may be stated in a more familiar form if 
the integrals are transformed by integration by parts 


(28) 


The results are then 


: 8 “es (x) dx 
C. = 24a + [ a es 
c? J—-c/2[1 — (2x/c)| V1 — (2x/c)? 


. T 8 {2 — n(x)xdx (29) 
C y= oa = = 
2 c8J/-c/2V 1 — (2x/c)? 


These are the well-known Munk integrals. 
The lift distribution can be obtained from Eq. (3). 
In a linearized form this can be written as 


p = —pU(u — U) (30) 


where the pressure is taken as zero far from the airfoil. 

If this result is combined with Eqs. (20) and (26), it is 

seen that the pressure distribution over the airfoil is 
ee 


es Co 6 . : 
—pU (a — tan sin 76 31) 
we. wae ; = 


This expression can be used to calculate hinge moments 
or, if the lift distribution is given as in the inverse prob 
lem, Eq. (31) determines the required camber line. 


given by 


»h = 
r= 


THE EFFECT OF THICKNESS 


The analysis of an airfoil of small but finite thickness, 
as in Fig. 1, can be carried out in a similar manner. 
For this purpose it is convenient to consider the airfoil 
as being built up with a thickness function 7, and a 
camber function n, where 


m = */2(m — ) 


'/o(m + 12) (32) 
With negative values of @ corresponding to the bottom 
side of the airfoil and positive values to the top, 0n/Ox 
may be expanded in the Fourier series 


and n, = 


re) Co vs ao : al 
T= + > cn cos nO + >) d, sin n0 (33) 
Ox 2 n=1 n=1 


where 


$ 
II 


2/4 f 2n./22) cos né dé (34) 
/0 
and 


d, = 2/x fn /ds) sin n6 dé 35) 
0 
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The flow around the airfoil is then 
: a 2U C 
U(1 — ta) — (-a+ ) - 
ae 2 


Ui Nat+ US dc- 


n=1 n=1 


on =< 
“e —_— 


(36) 


From this, it is apparent that the thickness function 
has no effect on the lift distribution or lift or moment. 
These are all determined by the camber function. The 
lift, moment and lift distribution are thus given by 
Eqs. (28) and (31) if they are applied to the camber 
function. 

The pressure distribution due to the 
function is found from the linearized pressure equa- 


thickness 


tion to be 


—pU? >> d, cos né (37) 


n=1 


p = 


This can be added to the pressure distribution due to the 
camber, Eq. (31), to obtain the resulting pressure dis- 
tribution. 


NUMERICAL EXAMPLES 
The results of the analysis will now be applied to 
several examples. 
Example 1 


The first problem to be investigated is the determina- 
tion of the lift, pitching moment and lift distribution of 
a* parabolic arc airfoil (Fig. 4) at zero angle of attack. 











Fic. 4. The parabolic arc airfoil. 


for this case 


n = fc{l — (2x/c)?| (38) 
ind 
On/Ox = — 2fc(2x/c)(2/c) = — 4f cos 6 (39) 
From Eqs. (24) and (27), 
c = — 4fand bd,” = 4f (40) 
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All other coefficients are zero. The lift and moment 


coefficients are 


C, _ Afar 
> (41) 
Cc u = O 
From Eq. (31), the pressure distribution is 
pb = (—pU”")(4f sin 6) (42) 


This shows that the lift is elliptically distributed across 
the chord. 


Example 2 


As a second example, a camber function will be calcu- 
lated which will produce a lift distribution across the 
chord which is constant for the leading half of the chord 
and is zero for the trailing half. For this lift distribu- 
tion the pressure on the rear half of the airfoil is zero, 
and on the top of the front half it is 


b = — Ci(p/2)U? (43) 


This pressure distribution must be an odd function of 
the angle 0, and, if it is expanded in a sine series as in 
Eq. (31), the coefficients are found to be 


L T 
Cc. = £ cos nr — COS N— >? (44) 
nr | 2) 
The slope of the camber function is then 
On/Ox = >> c, cos n (45) 
n=1 


If this is integrated term by term using the condition 
that » = 0 at the trailing edge, it is seen 


. 1 — cos (n + 1)0 
+¥a(! ~ 
X n+ 1 


\@ 
Zh ) | (46) 


In this case and in most cases where the pressure dis- 


1 — cos 2 


| 
. Aecaiilinas C1 
4 2 


1 — cos (n — 


n— | 


tribution is geometrically simple the series may be 
summed and the camber function may thus be expressed 
in a closed form. If the coefficients c, of Eq. (44) are 
inserted in Eq. (36), it is easily seen that 


; ux | | = 
w= U — Ui~"| '/2log{ 1+ — ) — log(1+- (47) 
7 ‘ye 3 
On the unit circle where ¢ = e”’, this becomes 

, ef cos 6 

w= U — Ui— log (48) 
2a 1 + cos 6 

and the slope of the camber function is 
On = Cr on cos 6 (49) 

Ox 2Qr 1 + cos 6 


This can be readily integrated and the camber function 
is found to be 
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—— 2 log 2 + cos 6 log jcos 6} — 

Cic : 


(50) 


(1 + cos @) log (1 + cos 6) 
This camber function is shown in Fig. 5. It may be 
noted that the s-curvature of the airfoil is just sufficient 
so that the pitching moment about the quarter chord 


point is zero. 
REFERENCES 


1 Birnbaum, W., Die Tragende Wirbelfladche als Hilfsmittel zur 
Behandlung des ebenen Problems der Tragfitigeltheorte, ZAMM, 
Vol. 3, pp. 290, 1923. 

2 Durand, W. F., Aerodynamic Theory, Vol. 2, Julius Springer, 
Berlin, 1934. 

3 Millikan, C. B., Am Extended Theory of Thin Airfoils and Its 
A pplication to the Biplane Problem, N.A.C.A. T.R. No. 362. 

* Keune, F., Die ebene Potentialstrémung um Allgemeine dicke 
Tragfltigelprofile, 1938 Yearbook of German Aero. Research, 
p.113. 


Letter to 


Dear Sir: 

In the paper by F. N. M. Brown, “‘The Space-Time Recorder” 
(June, 1942, JouRNAL), it is not stated what accuracy the photo- 
graphic type is capable of except to affirm its superiority to other 
types which are described as having an inertia error of the order 
of 0.1 to0.15g. One source of error not discussed by the author 
arises in taking off the displacements from the record. Using 0.1g 
as a desirable limit of error and !/2,9 of a second as the time inter 
val, we may calculate the displacement error which is 


6.1 xX 3x I. 


2 X 240 X 240 

This dimension, one-third of one-thousandth of an inch, is 
further reduced by the scale of the image. In the instrument de 
scribed by the author the reduction factor is one-tenth to one 
half. The probable. error in taking the displacements off the 
record is of the order of at least ten times this magnitude or a cor- 


responding error in the acceleration of 2g 
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(a) The lift distribution produced by the camber function 


(b). 


Fic. 5. 


5 Gebelein, H., Verallgemeinerung der Birnbaumschen Theorie 
ftir die Potentialstrémung um mdssig dicke Profile, 1938 Year 
book of German Aero. Research, p. I 27. 


the Editor 


One obvious cure for this error is to increase the time interval 
But this has its limitations since the acceleration calculated is t} 
average over the longer time interval. 

As a means of measuring airplane landing stresses, the space 
time recorder, along with the accelerometer, is of limited useful 
ness. Unfortunately, the airplane is not a rigid body and the 
simultaneous accelerations of its various masses vary widely 
Since the total landing force depends on the accelerations of all 
the masses, we cannot use a single acceleration for its calculation 
unless we know the shape to which the structure is distorting 
The accuracy with which we must know this distortion is in 
From it, 


cated by the preceding calculation. it is evident that tl 


assumption that equal accelerations occur at two points in 
structure is in error by 0.1g when a relative distortion at the 


points of 0.00033 in. occurs in !/249 of a second 


E. S. S 


JENKID 
Chief of Materials and Structures Resear 
Curtiss-Wrigl 


it Corporat 














A Simple Method of Applying the Compres- 
sibility Correction in the Determination 





of True Air Speed 


WM. C. SCHOOLFIELD* 
United Aircraft Corporation 


SUMMARY 


By reviewing the complete impact pressure equation, which 
describes the pressure exerted by an air stream upon coming to 
rest, it is shown that the approximating assumption that is in- 
volved in the customary method of determining the true air speed 
of an airplane at altitude is no longer satisfactory for high-speed 
The compressibility error that is a result of this as- 
sumption is isolated, and a method is described whereby the neces 


airplanes. 


sary correction may be applied by constructing additional alti- 
tude curves of indicated air speed (V+/c) on the conventional 
air-speed indicator calibration chart. The true air speed at any 
altitude is then determined by the usual method- the 
known value of indicated air speed and atmospheric pressure and 


from 


temperature. The error that results if the compressibility cor- 
rection is not applied is shown to be +11 m.p.h. for a 400-m.p.h 
airplane at 25,000 ft. altitude. The method here presented for 
applying the compressibility correction is less cumbersome than 
other suggested methods and has the additional advantage of 
retaining the use of the customary air-speed indicator calibration 
chart. Finally, the possible effect of compressibility on the posi- 
tion error of a pitot-static head installation is briefly analyzed 
and discussed, and the difficulties encountered in determining 
true atmospheric pressure and temperature in high-speed flight 
are mentioned. 


INTRODUCTION 


a THE TITLE of this paper should more properly 
be ‘‘A Simple Method of Using the Complete Im- 
pact Pressure Equation in the Determination of True 
Air Speed’”’ rather than its present title, for, actually, 
the procedure suggested in the following discussion is 
simply the correct application of a well-known and 
well-established law. The word correction arises from 
the fact that an approximating assumption has been in 
general use for so long a time that the error involved 
has been generaily overlooked. Two assumptions are 
frequently made in dealing with air-stream velocities 
and impact pressures. One of these assumptions is 
that the impact pressure of the air upon coming to rest 
is proportional to the square of the velocity multiplied 
by the air density, and the second is that for a constant 
value of the impact pressure the product of the square 
of the velocity times the air density is always constant 
and independent of the separate values of pressure 
and temperature. It will be shown in the following 
discussion that the first of the above assumptions is not 
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involved in the conventional method of determining 
true air speed but that the error that has been generally 
overlooked arises from the second of these assumptions. 
In order to appreciate satisfactorily the situation it is 
desirable to review the complete impact pressure equa- 
tion and to examine the conventional procedure for 
calibrating the air-speed system of an airplane. 


THE IMPACT PRESSURE EQUATION 


The usual air-speed indicating system in present-day 
airplanes consists of a pitot-static head located in the 
air stream and generally on the forward part of the 
wing or fuselage, connected by two tubes to an air- 
The 


cockpit instrument contains a sealed diaphragm, suit- 


speed indicating instrument in the pilot’s cockpit. 


ably connected to an indicator hand, and is operated by 
the pressure difference between the total, or stop, pres- 
sure of the air stream and the static pressure, these 
pressures being transmitted by the pitot-static 
head. 

The difference between the total pressure and the 
static pressure will be called the impact pressure, and 
it becomes apparent that the air-speed instrument is 
nothing more than a differential pressure gauge designed 
to be operated by the impact (differential) pressure of 
the air stream. For an incompressible fluid the im 


pact pressure is given by the expression 


P:—- p= pV?/2 = q (1) 
where 
p, = total, or stop, pressure 
p = static pressure, 
p = air mass density 
V = velocity of air stream 
g = usual symbol for p\V?/2 


The differential pressure expressed by the above equa 
tion is usually defined as the ‘‘dynamic pressure,’ and 
it is to be noted that the dynamic pressure so defined 
is the impact pressure of a fluid stream predicated on 
The 
air of the atmosphere is not incompressible, however, 


the assumption that the fluid is incompressible. 
and therefore Eq. (1) must be modified. As shown by 
Zahm,' the total pressure exerted by a moving air 
stream at the stagnation point of a body immersed in 


the air stream is 








THE 


V2 k/(k — 1) 
: | (2) 


where the new symbol, k, is the ratio of specific heat of 
air at constant pressure to its specific heat at constant 
This equation is derived from Bernoulli's 
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(k — 1) 


p= p[it kp 


volume. 
equation for compressible flow with adiabatic compres- 
sion, and the validity of the formula has been demon- 
strated by numerous experiments. For more readily 
comparing this equation with Eq. (1) it may be put in 
the following form, giving an expression for the impact 


(k = 1) V2 k/(k — 1) 
p—p= eli + | 1) (a) 


In order to show further the significance of this ex- 
5 
pression the term on the right may be expanded by the 
binomial theorem, at the same time substituting for 
l/s r . . o . 
(kp/p)’* the symbol JV, yielding the following: 


V 2 
pi — p = (pV?/2) [ + ie) + 


2-2) (V\! 
SS ee (2b) 


In Eq. (2b) the symbol lV, represents the velocity of 
sound in air and is related to atmospheric conditions 


pressure: 


as follows: 


velocity of sound 


V. = (kp/p)'” =a/T = 


where 7 = absolute temperature and a = constant. 
It will be noted that the velocity of sound in air varies 
only with the temperature and, if 7 is in degrees centi- 
grade (and for k = 1.400), the value of the constant a 
is 44.85, giving the velocity of sound in miles per hour. 

Eq. (2b) shows that for compressible fluids the actual 
impact pressure is always greater than the ‘dynamic 
pressure,’’ and Eq. (2a) shows that the relationship 
between the impact pressure, which will be designated 
q-, and the dynamic pressure, g, is uniquely determined 
by p, the static pressure. (The specific heat ratio, k, 
is assumed constant over the normal range of atmos- 
pheric temperatures.) 


AIR-SPEED INSTRUMENT CALIBRATIONS 


At this point it appears desirable to explain how the 
dial of an air-speed indicator instrument is designed or 
calibrated. It is to be made clear that this immediate 
discussion applies only to the instrument itself and not 
to the overall flight calibration of the air-speed system, 
which includes the pitot-static head. The conventional 
differential pressure type of air-speed instrument can 
be designed to indicate the true value of the air speed 
only under one set of atmospheric conditions as deter- 
mined by temperature and pressure independently. 
The Air Services of the United States Army and Navy 
in 1925 adopted the adiabatic formula of Eq. (2a) for 
the calibration of air-speed indicators? and selected 
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standard sea-level atmospheric conditions for the defi- 
nition of the temperature and pressure to be used in 
the formula. Consequently, for the marking of the 
instrument dial, Eq. (2a), after the substitution of 
standard sea level density and pressure, may be used 
to calculate the differential pressure corresponding to 
any value of the air speed, V. The value of the con 
stant, k, is taken as 1.40 for air. A table of these 
differential pressures associated with various values of 
V at sea level is given in Table IV of reference 2, and 
it is not necessary to reproduce this table here. It will 
be noted that the compressibility of the air 7s taken into 
account in the calibration of air-speed instruments; 
however, difficulties immediately arise in making the 
proper interpretation of an instrument reading when 
this reading is obtained under atmospheric conditions 
differing from the standard sea-level conditions for 
which the calibration is made. Of course, it is appreci 
ated that the air-speed instrument may be considered 
merely as a pressure gauge, and for a known value of 
the impact (differential) pressure as indicated by the 
pointer hand, Eq. (2a) may be used to solve for the 
true speed under any conditions of pressure and tem 
However, this process is cumbersome and 


perature. 
The method 


tedious and is to be avoided if possible. 
by which such calculations may be avoided will be 
shown later, but before continuing it is desired to de- 
fine the term true instrument air speed. ‘True in- 
strument air speed’ will be defined simply as the 
reading of a perfect instrument in a perfect installa- 
tion—that is, the reading of a correctly calibrated 
instrument when subjected to the true impact pressure 
of the air stream. This definition merely fixes the 
relationship between true instrument air speed and the 
impact pressure according to the standard formula Eq. 
(2a), which has previously been explained. 


OVERALL CALIBRATION OF AIR-SPEED SYSTEM 


The overall calibration of the air-speed system of an 
airplane includes the functioning of both the pitot- 
static head installation and the cockpit instrument. A 
conventional air-speed indicating system, in general, 
never gives the correct value of the indicated air speed, 
so an overall calibration of the system is always neces- 
sary. It will be immediately necessary to define the 
term “indicated air speed.’ The usual and simple 
definition of indicated air speed, which will be used 
here, is that indicated air speed is the product of the 
true air speed and the square root of the relative density 
ratio—that is, 
indicated air speed = V/p/p, = VVo (4) 
Rather than confuse the discussion with additional 
symbols, the symbol for indicated air speed will be its 
equivalent in symbol form—viz., Vo. As a matter 
of interest, the author has been unable to find a formal 
definition of indicated air speed in N.A.C.A. literature. 
The definition used above, however, is the same as that 
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given in recent C.A.A. publications, and it appears de- 
sirable that this definition should be standardized by 
all agencies concerned. It will be brought out later, 
however, that even a perfect air-speed indicator instru- 
ment in a perfect installation will not indicate true 
indicated air speed except under the special condition 
of sea-level pressure. 

If the above definition of indicated air speed is 
accepted, and the differences between instrument air 
speed and true indicated air speed are investigated, then 
the errors in a conventional installation may be classi- 
fied as follows: (1) position error of the pitot-static 
head, (2) calibration error of the air-speed instrument, 
(3) compressibility error, and (4) lag and acceleration 
error. It is necessary to understand these usual errors 
in order to fully appreciate the possible effects of com- 
pressibility. Item (4), however, which involves changes 
in altitude and changes in speed, will not be considered 
here. The present discussion will be confined to level 
flight, steady state conditions. 

(1) The position error, or installation error as it is 
sometimes called, is concerned with the pitot-static 
head and its location. The error arises from the fact 
that either the pitot tube may not be transmitting the 
true total pressure of the air stream or the static tube 
may not be transmitting the true static pressure of the 
air stream; or the error may be due to a combination of 
these conditions. A properly designed pitot-static 
head, when placed in a uniform air stream with the 
flow parallel to the axis of the tube, will normally trans- 
mit the total pressure and the static pressure of the air 
stream with satisfactory accuracy. And, furthermore, 
such a head is reasonably insensitive to small angles of 
pitch up to 8 or 10 deg. This is equivalent to saying 
that the constant of a well-designed head may be con- 
sidered unity over a limited range of pitch angles. A 
pitot-static head mounted on an airplane, however, is 
of necessity so close to either the wing or fuselage or 
some other part of the airplane that it is definitely in a 
region of flow which is influenced by the adjacent part 
of the airplane. The local flow at the pitot-static head, 
in general, does not have either the same velocity or the 
same static pressure as the free air stream, and the flow 
is usually at some angle to the axis of the tube. The 
pitot orifice, however, transmits the total pressure of 
the air stream, which is the sum of the impact and the 
static pressures; since this sum is always constant in 
streamline flow—i.e., flow for which Bernoulli’s equa- 
tion is applicable—the pitot tube gives a correct indica- 
tion of the total pressure regardless of the local value of 
the velocity and of the static pressure. The local static 
pressure that exists at, and is transmitted by, the static 
tube is, in general, not the true static pressure of the 
undisturbed air stream. Consequently, the differential 
pressure transmitted by the pitot-static head is usually 
not the true impact pressure of the free air stream. In 
the high-speed range of airplanes the local angle of at- 
tack is generally small and, if the tube is relatively in- 
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sensitive to pitch angles, no appreciable error will result 
from this source. Therefore, it may be concluded that 
the primary source, and usually the only source, of the 
installation error is the incorrect static pressure that is 
transmitted by the static tube from the pitot-static 
head. 

The above discussion must be considered to apply to 
the operation of the pitot-static head in its subcritical 
flow régime—that is, below the speed at which a com 
pressibility shock wave is formed near the nose of the 
tube. Experiments on typical pitot-static tubes** 
show that the correct total pressure is transmitted by 
the pitot orifice up to the sonic velocity, but the pres- 
sure at the static orifices begins to depart appreciably 
from the correct static pressure at a Mach number in 
the range from 0.72 to 0.80. Reference 4 satisfactorily 
correlates this error in static pressure at high Mach 
numbers with the formation of a shock wave near the 
nose of the tube and forward of the static orifices. It is 
felt that this error should be classed as a change in the 
constant of the pitot-static head rather than as an in- 
stallation error, and although it is definitely a compressi- 
bility phenomenon it should not be confused with the 
subcritical compressibility error, item (3), which will be 
explained later. At the present time the level flight, 
maximum speeds of airplanes are not high enough to 
exceed the critical speed of the pitot-static head; how- 
ever, the maximum diving speeds of modern fighter and 
pursuit airplanes are now reaching a range where pos- 
sible errors from this source must be given considera- 
tion. Since the discussion of this paper, as stated 
previously, is confined to level flight conditions, the 
constant of the pitot-static head will be taken as unity 
and no further reference will be made to errors arising 
from supercritical velocities. 

(2) The calibration error of the air-speed instrument 
is simply the error that may exist between the specified 
position of the indicating hand (according to the cali- 
bration formula) and the actual position of the hand for 
a given differential pressure on the instrument. This 
may be due to small errors in marking the dial, or it 
may be due to faulty mechanical operation of the in- 
strument mechanism. For the sake of completeness, 
temperature errors should also be mentioned here. 
Even though modern instruments are designed with 
temperature-compensating devices, there may be a 
small but appreciable change in calibration when the 
temperature of the instrument is varied over a wide 
range. 

(3) The compressibility error is defined as the dif- 
ference between the true instrument air speed and the 
indicated air speed, V+/o. Perhaps this difference 
should not be classed as an error at all, for it is evident 
that if the overall errors of an air-speed system had 
been defined as the difference between the instrument 
air speed—i.e., the observed reading of a given instru- 
ment—and the true instrument air speed this so-called 
compressibility error would not have appeared at all. 
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Nevertheless, since indicated air speed as herein defined 
is a simple quantity from which the true speed may be 
readily obtained (after the density ratio is known), the 
determination of this quantity will be considered the 
primary aim of this paper. Consequently, the re- 
mainder of this discussion will deal with the isolation 
of the compressibility error and a means of readily ap- 
plying it in the determination of indicated air speed. 


THE COMPRESSIBILITY ERROR 


Indicated air speed is a direct function of pV*/2 or q, 
since 
Vo = VV 0/m = V24/p0 (5) 
where po = standard sea-level density = constant. 
Now, referring to Eq. (2a), it will be recalled that when 
standard sea-level values of pressure and density are 
substituted in this equation it becomes the formula 
for the design calibration of air-speed instruments. 
Consequently, under sea-level conditions of pressure 
and temperature it is apparent that true air speed and 
true instrument air speed will be identical. Also, since 
at sea level o is equal to unity, it is evident that indi- 
cated air speed will have the same value as true air 
speed. Therefore, it develops that under standard sea- 
level conditions true air speed, true instrument air 
speed and indicated air speed are all the same. Now, 
let sea-level pressure be retained and investigate what 
changes take place when the atmospheric temperature 
The investigation will be made at a constant 
therefore, since the true 


is varied. 
value of the impact pressure; 
instrument reading is a function only of impact pres- 
sure, the true instrument air speed will remain constant. 
It is obvious that if the temperature is varied and the 
pressure held constant the density, p, must vary. Eq. 
(2a) shows, however, that if p and the impact pressure 
are to remain constant the product p\)* must remain 
constant. And since indicated air speed is a direct 
function of pV", then the indicated speed will remain 
constant also. It will be seen that the true speed as 
well as the density ratio must vary, however. This 
situation may be explained in a different way by the 
use of Eq. (2b). If the impact pressure—that is, 
bp. — p— and the quantity pV* both remain constant, 
then it is apparent that the term jn the brackets on the 
right must also remain constant and that the ratio of 
V/V, must remain unchanged. If pV? remains con- 
stant with varying temperature, then )* must vary as 
7’, the absolute temperature, since p is proportional to 
p/T. This variation, however, is exactly the variation 
shown for V’, by Eq. (3). Consequently, the ratio of 
V/V, in the bracket of Eq. (2b) is unchanged, which is 
in accord with the indications from Eq. (2a). The 
main point to be noted here is that at a constant value 
of the impact pressure, the Mach number, i.e., the 
ratio |/V., is independent of temperature; conse- 
quently, if the Mach number changes at all it must vary 
with the pressure only. Before passing to the case of 
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varying pressure, it will be observed that under the 
conditions of sea-level pressure the true instrument air 
speed and the indicated air speed are the same, regard 
less of the value of the temperature. 

For the case of changing static pressure, p, which 
represents a variation of altitude, it is desirable Lo ob- 
tain some convenient relation between true instrument 
air speed and indicated air speed. This relation will 
be sought at a constant value of the impact pressure, 
and it is to be noted that if the impact pressure is con 
stant the true instrument air speed will remain con 
stant. Reference will again be made to Eq. (2a). It 
has been shown that indicated air speed is a simple func- 
tion of pV*. Eq. (2a) may be rearranged for a solution 
of pV*, at the same time substituting for the impact 


pressure, p, — p, its symbol, g,, and the following is 
, kp | (2: 
V2/2 = = 
pve, i-1L + 


(k — 1)/k 
l ) — | | (2c 
and the indicated air speed is 


{2k (2 iin 
| 7” ee 
la@-itelw 


(2d) 


obtained: 


Vv -* = 


Let Vi~/ o; represent the indicated air speed under sea 
level conditions of pressure. Then 


2kpo 


Vivo = \(b — }) 


and the ratio of indicated air speed at any pressure, p 
to the indicated air speed at sea level pressure, po, is 


Age" 


Vo ] 
= ———— (f 


Me hee] 


where the exponent 0.2857 is the value of (k — 1)/k 
when & is equal to 1.40. The values of this equation 
have been computed for various values of g, and Pp, 
covering a range of g, represented by true instrument air 
speeds from 200 to 600 m.p.h. and a range of p repre 
sented by pressure altitudes from sea level to 40,000 ft. 
The calculated values are presented in Table 1 to four 


iO 


significant figures, which permits the calculation of the 
indicated air speeds at altitude to '/1) m.p.h. when so 
desired. Table 2 is an extension of the data of Table |, 
showing values of indicated air speed at altitude for 
constant values of the impact pressure, as represented 
by various indicated air speeds at sea level. It will be 
noted that blank spaces are left in both tables in the 
lower right-hand corner in the high-speed, high-altt 
tude region. This is done because the true air speeds 
represented by these spaces are above the speed of 
sound, in which region the applicability of the basic 
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TABLE 1 


RATIO OF INDICATED AIR SPEED AT ALTITUDE 


Standard 


Pressur¢ 


ro INDICATED AIR SPEED AT SEA LEVEL, 
INDICATED AIR SPEED = 


AT CONSTANT VALUES OF THE IMPACT PRESSURI 


Vy o 


Altitude Values of Indicated Air Speed at Sea Level. Vi~/o,;, M.P.H 
(Ft.) 200 250 300 350 400 450 500 550 600 
Sea level 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 
5,000 0.9982 0.9975 0.9963 0.9953 0.9938 0.9923 0.9907 0.9890 0.9873 
10,000 0.9960 0.9939 0.9916 0.9889 0.9859 0. 9827 0.9793 0.9758 0.9723 
15,000 0.9935 0.9902 0.9862 0.9817 0.9769 0.9717 0.9664 0.9608 0.9552 
20,000 0.9901 0.9851 0.9791 0.9727 0.9666 0.9583 0.9507 0.9430 
25,000 0.9860 0.9790 0.9708 0.9619 0.9524 0.9426 0.9326 
30,000 0.9809 0.9713 0.9606 0.9492 0.9365 0.9240 
35,000 0.97438 0.9617 0.9478 0.9329 0.9176 
10,000 0. 9660 (). 9499 0.9323 0.9140 
TABLE 2 
VALUES OF INDICATED AIR SPEED AT ALTITUDE, AT CONSTANT VALUES OF THE IMPACT PRESSURI 
Standard 
Pressure — Indicated Air Speed at Sea Level Vivio, = M.P.H 
Altitude (Ft.) 200 250 300 350 400 450 500 550 600 
Sea level 200 250 300 350 400 450 500 550 600 
5,000 199.6 249.4 298.9 348.3 397.5 446.5 495.4 544.0 592.4 
10,000 199.2 248.5 297.5 346.1 394.3 442.2 489.7 536.7 583 .4 
15,000 198.7 247.5 295.9 343.6 390.7 437.3 483 .2 528.4 573.1 
20,000 198.0 246.3 293 .7 340.4 386 .7 431.2 475.4 518.6 
25,000 197.2 244.7 291.2 336 .7 380.9 424.2 466.3 
30,000 196.3 242.8 288 .2 332.2 374.6 415.8 
35,000 194.9 240.4 284.3 326.5 367 .0 
40,000 193.2 237.5 279.7 319.9 
Eq. (2 is doubtful.’ : I he data of Table 2 are plotted 500 
in Fig. 1 as indicated air speed versus true instrument 
air speed, the latter being identical to Viv/o;. It will x 
be observed in Fig. 1 that true instrument air speed > 
has a different calibration—in terms of indicated air lo 
speed—for each value of the pressure altitude, which > 
represents a certain value of the static pressure, p. The O 400 
magnitude of the effect of pressure altitude on indicated 
air speed, for a given value of instrument air speed, will o SEA LEVEL ve 
be shown by an example. For the case of 400 m.p.h. a , 
true air speed at 25,000 ft. (¢ = 0.4480) the indicated = 
air speed is 268 m.p.h., and from Fig. 1 the true in- < f SOOO FT. 
strument air speed is 275 m.p.h. If, however, for an o 300 JY 25000 FT. 
instrument air speed of 275 m.p.h. the indicated air ra a 35000 FT. 
speed were read from the sea-level curve, giving an 5 | iy IOOOO FT ALL ALTITUDES 
indicated air speed also of 275 m.p.h., the true speed Pe - c 
at 25 io gles Anenitestines at ‘it at h eam 4 jy wee Mn Sa 
re ite: : OPN * WY 3SOOOOFT. ALTITUDES 
in an overestimation of the true speed by 11 m.p.h. ’ Cam 
This error, as indicated by the curves of Fig. 1, is less £200 AODOOFT. 
for lower altitudes and speeds but greater for higher 300 400 


altitudes and speeds. 

It will be observed that the calculations of Tables 1 
and 2 are based on a value of 1.40 for the specific heat 
ratio, k. The value of k, for dry air, found in various 


sources ranges from 1.40 to 1.41, and it appears that 
this ratio is actually a function of temperature, though 
the variation is extremely small (of the order of a few 





| TRUE INSTRUMENT AIR SPEED MPH. 


Fic. 1. Indicated air speed vs. true instrument air speed for 


various pressure altitudes. 


tenths of 1 per cent) over the range of temperatures 
encountered in the atmosphere. The value of 1.40 
used ‘here is consistent with that used in N.A.C.A. 
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publications. Check calculations for the case of 275 
m.p.h. true instrument air speed show that a change in 
k from 1.40 to 1.41 results in an increase in indicated air 


speed of less than '/1) m.p.h. at 25,000 ft. 


FLIGHT CALIBRATION OF AIR-SPEED SYSTEM 


It has previously been shown that the several errors 
that may exist in an air-speed system make it necessary 
to obtain an overall calibration that will include the 
errors of the pitot-static head as well as the errors of 
the cockpit instrument. Such calibrations are generally 
made in flight, and the usual procedure is to fly the 
airplane over a speed course at various constant in- 
strument air speeds and at the same time make observa- 
tions, through facilities that are independent of the 
airplane's conventional air-speed system, from which 
the airplane’s true air speed, as well as atmospheric 
pressure and temperature at the flight level, may be 
determined. From these data it is possible to compute 
true indicated air speed—that is, Vo. If desired, 
the true impact pressure of the free air stream may also 
be calculated by means of Eq. (2a). The usual pro- 
cedure, however, is to plot true indicated air speed 
against observed instrument air speed and obtain a 
curve that will give a relationship between these two 
speeds for any value of the instrument reading. A 
typical curve of this type for a monoplane on which the 
pitot-static head is mounted forward of the leading 
edge of the wing is shown in Fig. 2, labeled speed course 
calibration. The existence of an appreciable position 
error will be noted if it is assumed that the instrument 
error is small, as is usually the case. After such a curve 
of indicated air speed versus instrument air speed is 
obtained, it has been customary in the past to assume 
that this curve was applicable for flight at all altitudes 
and that for any condition of level flight the true air 
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Fic. 2. Typical calibration chart for an air-speed system; in- 
dicated air speed vs. observed instrument air speed. 
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speed could be determined simply on the basis of the 
density ratio and the indicated air speed from this curve. 
As soon as it is appreciated, however, that the cockpit 
instrument is operated essentially on the basis of an 
impact (differential) pressure and that a constant value 
of the observed instrument air speed indicates a con- 
stant value of the impact pressure of the free air stream, 
it becomes apparent from the previous discussion of the 
relationship between instrument air speed and indi 
cated air speed that the above assumption is in error. 
Although this error was negligible for slow-speed, low- 
altitude airplanes, it cannot be neglected in the ac- 
curate speed determination for high-speed airplanes 
operating at high altitudes. 

It remains now to show how the compressibility cor- 
rection to indicated air speed which has previously been 
computed may be applied to the speed course calibra 
tion of the overall air-speed system. The correction is 
made at a constant value of the instrument air speed, 
since this implies that the free stream impact pressure 
is constant. In order to apply the relations that have 
previously been computed it is necessary to know either 
the true instrument air speed or its equal, 11/«, 
since either of these quantities fixes the value of the 
If the speed course calibration is ob- 


impact pressure. 
pressure, then 


tained under conditions of sea-level 
Viv o, is immediately known, since the indicated air 
speed computed in the usual manner is identical to 
Viv/o, in this case. Therefore, the: data of Table 2 
may be used directly, and curves of indicated air speed 
at other pressure altitudes may be constructed on the 
usual calibration graph by plotting the points from 
Table 2 at the values of the abscissa (instrument air 
speed) which correspond to the tabular values of 11+ o; 
as indicated by the speed course calibration curve. 
Such a process has been carried out in Fig. 2 on the 
basis of the speed course calibration there shown, and 
calibration curves have been obtained for pressure 
altitudes up to 40,000 ft. in 10,000-ft. increments. 
These additional altitude curves of indicated air speed 
versus observed instrument air speed represent cali- 
bration curves that would have been obtained if a speed 
course calibration could have been made at each of 
these pressure altitudes and are to be used as such. 
Therefore, the procedure for determining the true air 
speed at any altitude for a given value of the observed 
instrument air speed is similar to that used in the past 
except that care must be taken to interpolate for the 
indicated air speed at the proper pressure altitude 
that is, at the pressure altitude at which the instrument 
air speed was observed. The true air speed is then 
readily computed after the density ratio is determined 
from observations of atmospheric temperature and pres- 
sure. 

In cases where the speed course calibration is ob- 
tained under atmospheric pressure conditions that vary 
appreciably from standard sea level, a correction may 
be necessary for the difference between Lhe speed course 
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value of Vo and the corresponding values of 1/0. 
This probably may be done most easily by cross-plotting 
the data of Table 1 to obtain the ratio of V+/o to 
Vivo, at any altitude. The process to be followed is 
essentially the same as that followed in constructing 
altitude curves from one obtained at sea level. It will 
be noted from an inspection of Table 1, however, that a 
relatively large departure from sea-level pressure will 
be required before this correction becomes large enough 
to be necessary. 

Overall calibrations of air-speed systems are some- 
times made by other means than flying the airplane over 
a speed course. One of these methods uses a special 
pitot-static head that is suspended at the end of a long 
cable at a sufficient distance below and behind the air- 
plane to be free from important interference effects. 
Since the primary cause of concern is usually only the 
static pressure, the suspended head may be only a static 
tube, in which case the total pressure is measured by a 
pitot tube directly attached to the airplane. In either 
instance it is presumed that the true impact pressure of 
the free air stream is correctly obtained. In this 
method of overall calibration, probably the simplest 
procedure is to determine true instrument air speed 
from the measured values of the impact pressure and to 
plot this speed against the observed reading of the 
normal cockpit instrument. Then, since true instru- 
ment air speed is identical to V1+/a;, the procedure for 
constructing curves of indicated air speed for various 
pressure altitudes becomes exactly the same as that 
described for the speed course calibration made at sea 
level. 

The author is familiar with two other methods® ® of 
handling the compressibility correction in determining 
true air speed; however, it is felt that the method here 
suggested is less cumbersome than other methods and 
has the additional advantage of retaining the customary 
air-speed calibration chart used in the past. Further- 
more, the magnitude of the correction involved is always 
readily apparent and the required interpolation is a 
minimum. 


EFFECT OF COMPRESSIBILITY ON POSITION ERROR 


In the foregoing discussion it has been assumed that 
for a given cockpit instrument air speed the true im- 
pact pressure of the undisturbed air stream was always 
constant, regardless of the altitude. This statement, 
of course, is restricted to a given airplane at a constant 
gross weight and assumes no change in the calibration 
of the air-speed instrument itself. In the general case 
where a position error of the pitot-static head exists, 
this assumption implies that the position error is con- 
stant as long as the instrument air speed is constant. 
Stated another way, the implication is that even though 
the differential pressure transmitted by the pitot-static 
head is not the same as the true impact pressure of the 
free stream these two differential pressures always bear 


a fixed relation to each other. The position error, as 
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previously explained, is generally due to the interference 
of the airplane wing or fuselage on the local static pres- 
sure at the openings of the static tube, and this inter- 
ference is a function of the lift coefficient or angle of 
attack and may also be a function of Mach number. 
The usual assumption in the past has been that as long 
as the instrument air speed remained constant the angle 
of attack and the lift coefficient, as well as the flow pat 

tern about the airplane, all remained unchanged. As 
a matter of interest, however, it can be shown that when 
the instrument air speed remains constant and the 
altitude is varied the lift coefficient, as well as the angle 
of attack, actually undergoes a slight change. This is 
due to the fact that the Mach number changes with 
altitude. Assuming for the present a constant value of 
the impact pressure and a change in standard altitude 
from sea level to 25,000 ft., an investigation reveals the 
following typical relations: 


Standard Altitude Sea Level 25,000 Ft 


True instrument air speed 275 m.p.h 275 m.p.h 
268 m.p.h 
400 m.p.h 


693 m.p.h 


Indicated air speed, Vo 275 m.p.h 
275 m.p.h 
761 m.p.h. 

0.361 0.577 


True air speed 
Velocity of sound 
Mach number, 


q-/ 1.033 1.086 
1/V/1 — M? 1.072 1.226 
Lift coefficient 0.1550 0.1635 
Slope of lift curve, dC, /da 0.0700 0.0800 
Angle of attack, absolute 2.21 deg. 2.04 deg. 


A wing loading of 30 Ibs. per sq.ft. is assumed, and the 
lift curve slope of 0.0700 is arbitrarily chosen for the 
sea-level case and then increased in proportion to the 
respective values of 1/+7/1 — M?*. It will be observed 
that the lift coefficient has increased 5.5 per cent, while 
the angle of attack has been reduced by 7.7 per cent; 
furthermore, these changes are accompanied, or caused, 
by a change in the flow pattern about the airplane. 
These changes are obviously small in absolute magni- 
tude; nevertheless, they do indicate that there may 
also be a change in local static pressure at the pitot- 
static tube relative to free stream static pressure. In 
order to form some idea of the magnitude of the possible 
change in local static pressure at the static tube openings, 
and the effect of this change on instrument and true air 
speed, it will be assumed that the static pressure error 
has been found at sea level to correspond to a Ap/g of 
+0.10, which is not an abnormal error for a pitot-static 
tube mounted ahead of a monoplane wing. The pres- 
sure error is here given in terms of the usual pressure 
coefficient as a function of g; however, it is well known 
that pressure coefficients of this type are a function of 
Mach number. On the basis of the present limited 
knowledge regarding compressibilty effects, it is gen- 
erally assumed that for an aerodynamic body at a 
constant attitude relative to an air stream the stagna- 
tion pressure coefficients vary with Mach number as 
q-/q, While all other surface pressure coefficients vary 


approximately as | /1 — M?, where M is Mach 
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number. The impact pressure acting on the open end 
of the pitot tube is a stagnation pressure and is identical 
with g,. This point is the basis of all the previous dis- 
cussion of the compressibility error and is satisfactorily 
proved by experiment. Therefore, it becomes evident 
that if the static pressure error varies in any other 
manner except directly with g, the position error will 
change with altitude. Continuing the comparison be- 
tween sea level and 25,000 ft., if the static pressure error 
is proportional to g,, the error expressed as a pressure 
coefficient will become (1.086/1.033) (0.10) = 0.1051. 
In this case the position error remains constant. If, 
however, the pressure coefficient varies as the factor 
L, V/1 — M2, then the coefficient at 25,000 ft. becomes 
(1.226/1.072) (0.10) = 0.1144. The ratio of the coef- 
ficient in the latter case to that in the former is 
0.1144/0.1051 = 1.09. This means that the error in 
static pressure at the openings of the static tube has 
increased by 9 per cent, while the impact pressure on the 
pitot tube remained constant, which reduces the differ- 
ential pressure acting on the air-speed meter and con- 
sequently results in a reduction of instrument air speed. 
The reduction in differential pressure would amount to 
1 per cent, corresponding to approximately '/2 per cent 
reduction in instrument air speed, and this change would 
take place, according to the conditions of the investi- 
gation, at a constant value of the true impact pressure 
of the free air stream. Consequently, the true air speed 
would be underestimated by approximately '/2 per 
cent, and, although this is a small percentage, it amounts 
to 2 m.p.h. true air speed in this selected case. If the 
assumed error in the static pressure at the static tube 
had been negative, then the true air speed would be 
overestimated by 2 m.p.h. 

The above computations are obviously approximate 
and are intended only to point out the possibility of the 
small error that may be involved in the assumption that 
the true impact pressure is constant at all altitudes as 
long as the instrument air speed is constant, even though 
the pitot-static head has a position error. Under the 
circumstances, however, it appears that sufficient infor- 
nation to warrant any other assumption is not avail- 
able. Obviously, the smaller the position error, the 
smaller will be the error resulting from changing alti- 
tude and, if there is no position error, the point of this 


discussion disappears entirely. 
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DETERMINATION OF DENSITY RATIO 


Only brief mention has been made in this paper re 
garding the determination of the density ratio, ¢, which 
is required in computing the true air speed after the 
indicated air speed is known. The determination of 
the density ratio requires a knowledge of the true at 
mospheric pressure and temperature at the flight alti 
tude and, while this subject is not strictly within the 
scope of this paper, it is desired to point out that these 
two quantities are not easily obtainable with satis 
factory accuracy in high-speed flight. It is customary 
to connect the altimeter to the static line from the 
pitot-static head, and when this is done it immediately 
becomes apparent that the pressure altitude that is 
indicated by the altimeter is subject to the same pres 
sure error that is reflected as a position error of the 
pitot-static head. While the error in pressure altitude 
from this source is relatively unimportant with regard 
to determining the compressibility correction to indi 
cated air speed, this error can be of definite significance 
in the determination of the true density ratio. Also, 
the usual thermometer that is carried on an airplane for 
measuring the free air temperature is always subject to a 
temperature rise due to compressibility effects. And 
this temperature rise is generally large at high velocities. 
Unless these two usual errors are properly appreciated 
and allowed for, they can readily result in errors in the 
true air speed of the same order of magnitude as thos¢ 
resulting from the neglect of the compressibility error. 


REFERENCES 


1 Zham, A. F., Pressure of Air on Coming to Rest from Variou 
Speeds, N.A.C.A. T.R. No. 247, 1926. 

2 Beij, K. Hilding, Aircraft Speed Instruments, 
T.R. No. 420, 1932. 

3 Lock, C. N. H., and Hilton, W. F., Calibration of Standard 
Pitot-Static Heads in the High-Speed Tunnel, British Air Ministry 
R. & M. No. 1752, August, 1936. 

* Walchner, O., The Effect of Compressibility on the Pres 
Reading of a Prandtl Pitot Tube at Subsonic Flow Velo 
N.A.C.A. T.M. No. 917, November, 1939. 

5 Thompson, F. L., The Measurement of Air Speed of Airplane 
N.A.C.A. T.N. No. 616, October, 1937. 

6 Kotcher, Ezra, The Compressibility Factor in Converting Ir 
cated to True Air Speed, U.S. Air Corps Technical Report N« 
4515, February 29, 1940, unrestricted. 


N.A.C.A 














Propeller Forces Due to Yaw and Their 
Effect on Airplane Stability 
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ABSTRACT 


Forces produced by a propeller operating at an angle to the air 
stream or being subjected to certain displacement velocities with 
respect to the air stream were calculated by Harris and Glauert 
in 1918. Recent tests on a propeller of modern plan form has 
shown that the theory of Harris and Glauert when used without 
their empirical constant predicted only half of the force in yaw 
as indicated by the tests. A new theory is presented here which 
has improved the agreement between the predicted and measured 
values of one of the force derivatives for the propeller tested and 
is believed to be more desirable for design calculations. 

This method is based upon the same assumptions used by 
Harris and Glauert but includes the modified lift produced by an 
The equations 
they 
depend mainly upon the propeller blade form and clearly indi- 


airfoil when operating in an oscillating air stream 
are in a different form than given by Harris and Glauert 


cate that propellers having a wide root chord section will produce 
a greater side force due to yaw than would be realized with a pro- 
peller blade having a narrow root chord. This is an interesting 
conclusion, since a high-altitude pusher-type airplane that will 
probably require a propeller with a wide root chord will be 
favored from the stability point of view in comparison with a 
similar tractor-type airplane. 


INTRODUCTION 


if be EXISTENCE OF FORCES produced by a propeller 
when the thrust axis is at an angle with respect to 
the air stream or when it is being subjected to an angular 
velocity with respect to the air stream has been known 
for some time. The basic theory of this subject was 
developed by Harris! for the case of the propeller side 
force due to angle of yaw, and it was extended to include 
the effect of angular velocities by Glauert.”, This work 
was done in 1917 and 1918 and until recently has not 
been extensively used, since the forces in question have 
generally been of secondary importance. Present de- 
signs require more engine power, necessitating more and 
larger propellers that may be disposed in new positions 
on the airplane. These tendencies make the problem 
of airplane stability more dependent upon certain forces 
produced by the propeller. 

Recent wind-tunnel tests of a powered wind-tunnel 
model provided a means of measuring the side force and 
yawing moment produced by a propeller. A comparison 
of these measurements with the theory of Harris and 
Glauert showed that at some values of |’/nD the test 
results were approximately twice those indicated by the 
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theory. For this reason a new method of analysis has 
been used and new equations found for two of the more 
important force derivatives. The new method of analy- 
sis gives, for the side force due to yaw, better agree- 
ment with test results than the equation of Harris. 
As the new equations depend mainly on the propeller 
planform, they are believed more practical from the 
standpoint of design usage. 


EQUATIONS OF HARRIS AND GLAUERT 


As the original works of Harris and Glauert are not 
readily available, the basic equations of reference 2 will 
be summarized and presented in a form that is believed 
to be more convenient for airplane stability calculations. 
Glauert shows that the following forces will be produced 
by a propeller undergoing a displacement or motion of 
the thrust axis: (a) side force (normal to thrust axis) 
due to angle of yaw, (b) pitching moment due to angle 
of yaw, (c) side force due to an angular velocity of 
pitch, (d) pitching moment due to an angular velocity 
of pitch, (e) force along the thrust axis due to a change 
in the velocity in the direction of the thrust axis, (f) 
torque produced by the propeller due to a change in the 
velocity in the direction of the thrust axis. 

These have been rewritten and presented in Table 1. 
All of the force derivatives have been expressed in 
terms of airplane dimensions using N.A.C.A. nomen- 
clature and are written for a right-hand propeller having 


TABLE 1 


Propeller Force Derivatives Based on Theory of Harris and 


Glauert for Constant Speed Right-Hand Propeller 


Aerodynamic 
‘alue of R 


a 


Derivative 


1Cy D? f| Cp 1 /dCp 
(a) a = R, = ( : ) 
dy 3 ar\ J 2\dJ/s 
dCy D3 11 GC: 1 fdCy 
b) = R, — ( ) 
dy St ri J 2\dJ/B 


Equation 


dCy D3 (tt) 
Cc) = 4 R 
d(qb/2V) Sb rJ\dJ/8 
dCm Ds f, (dCy 
(d) - = R, : 
d(qb/2V) Sht 2J\ dJ /B 
iC, D? -4[Cr 1/dC 
(e) ; . wa = 2 R; = ( J ) : 
d(u/t ) s yf P Zt Ge Jt P 
dC; 
(f) = = zero 
d(u/V) 
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the airplane axes coincident with the propeller hub in 
order to simplify Table 1. In applying these formulas 
it is a simple matter to include the moments produced 
by the forces when the airplane axes are referred to the 
airplane c.g. In the case of a left-hand propeller Eqs. 
(b) and (c) will change sign, hence must be considered 
zero for the case of a dual rotating propeller. The cal- 
culations and formulas of reference 2 were, of course, 
for a fixed-pitch propeller but, with two exceptions, 
are applicable to a constant-speed propeller. For con- 
stant-speed propellers Eq. (f) becomes zero and Eq. 
(e), as given in Table 1, has been recalculated. All 
other derivatives may be applied either to a fixed- 
pitch or constant-speed propeller. In the equations of 
Table | for Ri, Re, etc., the subscript 6 or Cp on a deriva- 
tive indicates the variable held constant during the 
differentiation. 

By letting g, and ¢, represent a factor proportional to 
the torque and thrust per unit length along the pro- 
peller blade (Eqs. (13a) and (24a)) values of f and f; 
used in R, and R, are given by the integrals 


f = (2/1) fo'q,(R/r)*d(r/R)/ So'gd(r/R) (1) 
fi = So't(r/R)*d(r/R)/ Jo't,d(r/R) (2) 


The equations of Table 1, together with a propeller 
chart, will permit the numerical values of the propeller 
force derivatives to be readily calculated by graphically 
measuring the necessary propeller coefficient deriva- 
tives. The formulas of this table may be used, of course, 
to calculate the force derivatives about other airplane 
axes by a proper interchange of the axes. 

In this paper the force derivatives (a) and (d) will 
be discussed in detail and new values of R; and R, will 
be determined by a different method. In addition, 
experimental values of (a) will be compared with values 
computed by the Glauert formula of Table | and by the 
new formula. 


SIDE FORCE PRODUCED BY A PROPELLER IN YAW 


In this analysis it will be assumed that a right- 
hand propeller is being sideslipped to the right, giving 
rise to a negative angle of yaw of the thrust axis but 
having a positive sideslip velocity. The rear view of 
this propeller is shown in Fig. 1.. The vector marked 
— VW indicates that the propeller hub is moving to the 
right normal to the airplane flight velocity at a side 
velocity — VW where W is the angle of yaw measured 
in radians. It will be assumed that all radial velocity 
components at the propeller blade may be neglected, 
which has been the assumption of all previous analyses. 
Fig. 2 shows the velocity vectors acting on the blade 
element indicated in Fig. 1 for the propeller blade in 


positions A and B. For the blade in position A no 


effect of the side velocity is experienced; therefore, the 
force coefficients will be denoted by the subscript » and 
will correspond to those for zero yaw. 

During the motion from A toward B it will be as- 
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sumed that the lift coefficient changes in magnitude and 
direction but that the drag coefficient changes only in 
its direction. 
zero yaw condition is unchanged. 
be the changes in the force coefficients of the blade ele 
ment normal to the radius in the plane of the propeller 
and normal to the plane of the propeller, respectively. 


Then by taking Aa negative as shown in Fig. 2 and 
of a small magnitude, cos Aa = 1, sin Aa = Aa, 


giving 


AC, =C, -— C,, 
= (C,, + mAa) sin (¢— Aa) + Cp, cos (6 — Aa) 


= (—C,, + m sin ¢) Aa (3 
Cz =? Cr 
(CL, + mAa) cos (¢ — Aa) — 


AC, = 


where 


Before proceeding with the integration of these force 
coefficients over the blade, the value of m will be estab- 
lished based upon Theodorsen’s Eq. (18) of reference ». 





2mrn 


SECTION LOOKING FROM HUB 
TOWARD TIP FOR BLADE AT @=90° 


(C,, + m cos ¢) Aa (4 


Fic. 2. 


This implies that the profile drag for 
Let AC, and AC, 


(C,, sin @ + Cp, cos ¢) 


) 


Cp, sin (¢ — Aa)—(C;,, cos ¢— Cp, sin ) 


m = dC,/da 














PROPELLER FORCES 


A particular blade element in moving from A to B may 
be considered as an airfoil describing an oscillation 
normal and parallel to the direction of the undisturbed 
air velocity. Considering the oscillation normal to 2, 
the vertical displacement velocity v, of Fig. 2 may be 
written as 
v, = vW¥ sin’ ¢ sin 0 

Reversing the sign of Eq. (18) of reference 3 in ac- 

cordance with the present notation and neglecting all 


terms not involving / or its derivatives, the lift coef- 
ficient AC, becomes: 


AC, = —P/(p/2)b’v? = [2xF sin 0 + 
(rk + 2G) cos 6|W sin? @ (5) 

where 

P = lift forces given by Eq. (18) of reference 3 

b’ = chord of blade element 

k = pb'/2v 

F,G = circulation functions of 1/k (see references 3 

and 4) 


In developing Eq. (5) the circular frequency is p = 
2xn where » = propeller speed (r.p.s.), and pt = @, 
making the amplitude of the relative displacement 


hy = (b'’/2k)WV sin? o 


It will be shown below that Eq. (5) must be multi- 
plied by sin 6, and all terms having the product sin @ 
cos 6 will be zero when the average side force during 
one revolution is calculated; hence, at this point the 
cos 6 term will be dropped. From Fig. 2 


Aa = YW sin? ¢ sin 0 (6) 
v = V/sin @ (7) 


But mAa of Egs. (3) and (4) is the lift coefficient 
AC, of Eq. (5); hence 


m = 2nrF (8) 


The value of F is a function of the wing aspect ratio 
and of 1/k = J/x(b’/D) sin ¢, where (b’/D) is the ratio 
of the propeller total blade element width b’ to propeller 
diameter D. In Fig. 3, curves of F are plotted for wing 
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aspect ratios of ~, 6, and 3, as taken from reference 4, 
and it is seen that for small values of 1/k the aspect 
ratio has little influence on F. In the present analysis 
the effect of aspect ratio will be neglected, all calcula- 
tions being based upon the curve for aspect ratio of 
infinity. 

The value of m given by Eq. (5) is that experienced 
by an airfoil describing an oscillation normal to v, but 
the blade element is also moving fore and aft, giving 
rise to a change in the of v of the 
amount 


magnitude 


Av = vW sin 6 sin ¢ cos ¢ 
or, combining with Eq. (6), 


Av/v = Aa/tan @ 


As this motion is in phase with the vertical oscillation, 
the value of m will be increased to include the effects of 
the change in the velocity v, which accompanies all 
changesin a. The value of m, however, will be referred 
to values of C, based upon v. If C,, is the lift coeffi- 
cient including the effects of the velocity v + Av but 
based upon v, then, as C, is proportional to v? and a, 
the following differential expression may be written 


AC,,/Cry, = (2Av/v) + (Aa/a) = 


Aa|(1/a) + (2/tan ¢)] (9) 
However, C,, = ma; hence, the modified slope of the 
life curve m, becomes 


m, = AC,,/Aa = m(1 + 2a/tan ¢) (10) 


For the purpose of this correction it is reasonable to 
assume that 


C,, = Cy, cos 
which permits a to be written in terms of C,, giving 
m, = 2nF(1 + 2C,,/2xF sin ¢) 
Substituting m, for m in Eq. (1) 
AC, = (C,, sin? @ + 2xF sin*®* ¢)W sin 6 


Referring to Fig. 1, the side force acting in the direction 
of the Y axis for the blade in position @ is 


VY = Jo*(b'/D)DAC,(p/2)v? sin 6 dr 
= (p/2)V*(D?/2)W sin 26 f,1(b’/D) X 


(C,, +22F sin ) d(r/R) (11) 


The above side force is periodic for a single blade, but, 
as the average side force is of interest here, the sin? 6 
term may be replaced by its average value of one-half 
as 


(: om) Jy” sin? 6d@ = ! 2 


When this is done the average side force coefficient 
Cy, based for convenience on the airplane wing area S, 
may be written for a propeller having B blades, as 


(BY)/(p/2)V2S = R,(D*/S)¥ 


Cy = 


(12) 
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where In order to get the first integral of R, in terms of the 


RK, = (8B 1) foi(d’ D)C,,d(r/R) + (x 2)B Sy '(b’ Dy xX 


F sin dd(r/R) (13) 


propeller thrust coefficient as measured for zero yaw, 
the thrust loading along the blade must be known if 


t, = (6'/D)C,,/sin? @ = factor proportional to the 





thrust loading (13a 
then the propeller thrust 7 is 
T = (p/2)V?(D?/2)BSfy't,d(r/R) = pn?D4C, 14 
giving 


(B 1) fo 't.d(r R) = C,/J 


+ 


The first integral of R; may now be written as 





(B 4) fo'(b’/D)C,,d(r/R) = 
B/4) Jo't, sin? ¢d(r/R) = M,C, 











° 2 4 6 8 
r/R 63 where 





| | | My, = (1/F*) SoM te/trm) sin? od(r/R) 
es mae | T | t T T 1 T So'(te/trm)d (r/R) 15 
The second integral of R, may be written as 
(7/2)BSy'(b'/D)F sin ¢d(r/R) = MoB 
where 


Me. = (b'/D)o-n(x/2 Si (b’/D)/(6’/D) 75r| X 


F sin od(r/R 16 


In order to determine .)/; and .Ws, the thrust distribu 





tion and blade plan form for the propeller in question 
must be known. For values of J in the flight range the 
term .J/,C; is generally small ‘n comparison to MW2B. 
This is fortunate since it enables an approximate thrust 
distribution to be used in determining .W, and then 
MJ may easily be calculated from the blade plan form. 
In order to present values of 1/;, three thrust dis- 
tribution curves (later also to be used as torque dis- 
oop} —_+_}__| __| | Ma | tribution curves) have been arbitrarily selected and 
| | | denoted by A, B, and C in Fig. 4a. These curves are 

| 








not associated with any particular propellers, with the 
exception that curve A represents approximately the 
thrust distribution for the propeller used in the model 
tests. Values of 1, for the thrust distributions A, B, 
and C are shown in Fig. 5. Likewise, .J/2 as plotted in 
Fig. 5 has been calculated for three arbitrary propeller 
plan form distributions denoted by I, II, and III of 
Fig. 4b. These curves are not related to the curves 
A, B, and C and are based upon (6/D)o 752 = 0.0817. 
Curve I represents an approximation to the plan form 
of the propeller tested. 

The final equation for the side force coefficient, con- 
sidered acting normal to the propeller thrust axis and 
based upon the airplane wing area, is 














Cy = R,(D? S)v 17) § 

















; 2 3 4 where 


R, = MC, + MB 1S) V 











PROPELLER 


PITCHING MOMENT DUE TO AN ANGULAR VELOCITY 
OF PITCH 


This aerodynamic derivative is calculated because 
it is useful for longitudinal stability calculations. It is 
obtained in a manner similar to the previous derivation, 
except that an axial change of velocity is analyzed in- 
When the propeller blade 
axis located in the 


stead of a transverse one. 
element is above the transverse Y 
plane of the propeller at the propeller hub, a positive 
angular-pitch velocity g (rad. per sec.) will cause the 
forward velocity V at the blade element to be reduced; 
for the blade positions below the Y axis an increase in 
V will be realized. Combining velocity vectors similar 
to those shown in Fig. 2, a blade element will experi- 
ence the following changes: 


Aa = (—q,/V) sin @ sin ¢ cos @ 
Av = q, sin 6 sin @ (19) 
Av/v = — Aatan 


In this case, when the blade element encounters a 
positive change in Aa, the value of Av is decreased; 
hence, the value of m must be reduced in order to ac- 
count for this variation in velocity. As an approxima- 
tion let 


C,, = C,,/tan @ = C,, cos ¢ 


C, = ma 
then by using the first expression of Eq. (9) 

n, = n[1 + 2(Av/v)(a/ Aa)} 

= 2rF[1 — 2C,,/2rF cos | 
Replacing this m, for m in Eq. (4) in order to account 

for the effect of the change in v, the force coefficient of 
the blade element force parallel to the thrust axis now 
becomes 
(21) 


AC, = (C,, — 24F cos ¢)(qr/V) sin ¢ cos ¢ sin 6 


The pitching moment for a propeller having B blades 
then becomes 

M = (p/2)(D®B/4) Jo'v2(b'/D) AC,(r/R) sin 6d(r/R) 

Substituting Eq. (21) in Eq. (22) and replacing sin* 6 
by its average value of '/2, the pitching moment coeffi- 
cient based upon the airplane wing area S and mean 
aerodynamic chord ¢ is 
Cn = M/(p/2)V?St 

= (D*/Sbt)(gb/2V) Jo'(b’/D) (r/R)? X 
(C,, — 2xF cos ¢)(1/tan @)d(r/R) (23) 


where 0 is the airplane wing span. However, tan ¢ 
may be expressed as 
tan @ = J/x(r/R) 
giving 
Cn = Ry(D4/Sbt)(qb/2V) 


where 


FORCES DUE TO 


YAW 169 


R, = (4b/8J) Jy'(b'/D)(r/R)®C,d(r/R) — 
rB/4J Jy'(b'/D)(r/R)*F cos ¢d(r/R) (24) 
The first integral of R, will depend upon the torque 
distribution along the propeller blade. 
Let 
sin* ¢) = factor proportional to 
(24a) 


gr = (6'/D)(r/R)(C 


the torque distribution 


then the torque Q on the propeller at zero yaw may be 
written as 
QO = (p/2)(V2D°B/4) f'¢,d(r/R) = pn2D*(Cp/2r) (25) 
giving 

So'qgd(r/R) = (4/2B)(Cp/J) 


The first integral of Ry then becomes 


(rB/SJ) Jy (r/R)*q, sin? dd(r/R) = NiCp 
where 
M = CJ So (r/R)2(G:/Grm) sin? ¢d(r/R) (96) 
So (Gr/Grm)d(r/R) 
The second integral of R, is written as 
—(x°B/4J) Jy (b'/D)(r/R)*F cos ¢d(r/R) —N.B 
where 
Nz = (m?/4J)(b'/D)o.23JSo ((b’/D)/(b'/D)o.7] X 
(r/R)*F cos ¢d(r/R) (27) 


The final equation for the pitching moment coeffi 
cient produced by a propeller due to an angular veloc- 
ity g about the transverse Y axis through the propeller 
hub is 
(28) 


Cy, = R,(D*/Stb)(qb/2V) 


where 


R, _ NiCp —- N.B (29) 


This expression has a form similar to that for the side 
force due to yaw. Using curves A, B, and C of Fig. 4a 
to represent the torque distribution along the propeller 
blade, values of N, have been computed and are plotted 
in Fig. 5. Values of Ne were likewise computed and 
plotted using the plan forms I, II, and III of Fig. 
4b. 


TEST RESULTS AND COMPARISON WITH THEORY 


The side force produced by a six-blade counter- 
rotating propeller has been determined from wind- 
tunnel tests on a powered wind-tunnel model. The 
effect of the running propeller was obtained by testing 
the model with and without the propeller. The model 
tested was such that any interference produced by the 
wing or fuselage was small. The side force acting at 
right angles to the propeller thrust axis was calculated 
and plotted against the angle of yaw, from which the 
values of dCy/dV were obtained from measurements of 
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M, M2 
Type B Type I 
B° J Cr Propeller Propeller 
30 1.33 0.075 0.13 0.040 
30 0.965 0.225 0.18 0.034 
38 1.5 0.192 0.137 0.042 
38 1.61 0.145 0.13 0.043 
38 1.04 0.315 0.172 0.036 
46 2.12 0.20 0.104 0.0475 
38 1.94 —0.035 0.113 0.046 
38 1.95 —0.035 0.113 0.046 
46 2.66 —0.10 0.082 0.051 
46 2.68 051 


0.082 0. 
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the average slope of these curves between 0 and 10 deg. 
angle of yaw. 

The propeller tested had a diameter of D = 2.25 ft. 
and a plan form similar to curve I of Fig. 4b with 
(6/D)o 22 = 0.0817. These tests results are tabulated 
in Table 2 in terms of R; along with the corresponding 
theoretic values of R;. 

The theoretic value of R, based upon Glauert’s 
equation (Eq. (a), Table 1), was obtained by using the 
propeller coefficients measured in the wind tunnel. 
A value of f = 1.44 was determined by using curve B of 
Fig. 4a to represent the torque distribution. Strip 
calculations have shown curve B to be a good repre- 
sentation of the torque distribution for this propeller. 

Theoretic values of R,; determined by the present 
theory are also given in Table 2. These values were 
obtained by using M/, and M, from Fig. 5 along with 
the measured thrust coefficients. 

A comparison of these results is shown in Fig. 6, 
which represents the ratios of the test values of R; to 
the values calculated by Glauert’s formula and by the 
present theory. 

It appears that an empirical correction factor of 
approximately 1.2 would bring the theoretic values into 
good agreement with the test values. Pending the 
accumulation of more test results, no final value of such 
an empirical factor can be recommended. 


DISCUSSION 


In this paper the theory of Glauert and Harris, 
which is based upon propeller coefficients, is compared 
with the present theory, which incorporates a correction 
to the blade element lift depending upon the theory of 
the oscillating airfoil. The present theory seems to be 
an improvement over that developed by Harris and 
Glauert and probably is a step toward obtaining a 
strictly rational formula. Before any conclusions on 
this point may be reached a further comparison with 
tests must be made. 

The present theory has the advantage that it de- 















































Glauert Present Glauert Present 
Test theory theory theory theory 
0.296 0.171 0.251 1.73 1.18 
0.296 0.191 0.244 1.55 1.21 
0.383 0.225 0.278 1.70 1.38 
0.395 0.221 0.277 1.78 1.42 
0.309 0.206 0.270 1.50 1.14 
0.383 0.268 0.306 1.43 1.25 
0.321 0.206 0.272 1.56 1.18 
0.309 0.206 0.272 1.50 1.14 
0.321 0.353 0.298 0.91 1.08 
0.358 0.353 0.298 1.01 1.20 
OPRESENT| THEORY | 
= I 
xGLAUERT!] THEORY f=1.44 
ore wr a 
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2 | x 4% | | 
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or ° 
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pends mainly on the propeller-blade plan form. As 
far as the side force due to yaw is concerned, it clearly 
indicates that propeller blades having large chords 
near the root will give the greatest side force due to 
yaw. The curves of Fig. 5 consist of values of M,, M2, 
N; and N, based upon three propeller-blade plan forms 
and thrust and torque distributions as shown in Fig. 4. 
The curves have been arbitrarily selected in order that, 
for any given propeller design, proper values of M/,, Mz, 
N, and N2 may be estimated permitting the side force 
due to yaw and pitching moment due to pitch to be 
readily calculated. 
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A Practical Method of Allowance for 
Shear Lag’ 


JOHN DE S. COUTINHO7 


Grumman Aircraft Engineering Corporation 


SUMMARY 

The calculation of stresses in a semimonocoque cantilever 
wing is discussed from an engineering-office point of view. 

The influence of various factors on the stress distribution in 
thin-walled structures is considered, especially the effects of 
shear lag for which a simple empirical method of allowance is 
suggested. The use of the method is illustrated in the example 
of a box-beam wing, and the results of the calculations are 
measurements on a full-size test speci- 


checked by stress 


men. 


INTRODUCTION 


I PRACTICAL AIRCRAFT DESIGN the strength-weight 
ratio is a factor of primary importance. Efficient 
design requires a small but reliable margin of safety in 
all elements. Hence, a reasonably accurate knowledge 
of the stress distribution is essential. 

In the commercial design of aircraft, however, the 
time element and the available mathematical skill 
often prohibit the preparation of intricate exact calcula- 
tions. The value of a method that is to be of practical 
use will therefore depend greatly on its simplicity, 
flexibility, and the time and labor involved in obtaining 
reasonable answers. 

This paper represents, therefore, an effort to con- 
tribute to the simple and practical empirical methods 
for the use of the skilled and experienced designer. 

IN A CANTILEVER BEAM 


THE STRESS DISTRIBUTION 


The design of a cantilever beam is based on the 
knowledge of the stress distribution in the various cross 
sections of the beam. 

Consider the stress distribution in a 
prismatic beam of narrow rectangular cross section, 
For a < b, the exact stress distribution is given 
normal stress, ¢ = My/I; shear stress, tT = 


cantilever 


Fig. 1. 
as follows: 
PQ bl. 

The normal stress distribution given in this ex- 
pression is plane and represents the simplest bending 
Stress distribution conceivable. The formulas are 
easily applicable, giving relatively simple calculations 
and quick results. The conception of planes before 
bending remaining planes after bending yields a de- 
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scriptive mental picture that is easy to visualize and 
to work with. 

For even slight deviations from the narrow rec- 
tangular cross section, however, the formulas, in all 
cases for which an exact solution is known, become 
considerably more complicated, though the actual dis- 
tribution may be similar. Since the possible gain in 
accuracy is out of proportion to the additional labor 
involved, the general practice is to assume plane stress 
distribution for all cases of bending. in the analysis of 
complex thin-walled structures this method may result 
in considerable error. The simple formulas are then 
applied to a modified and idealized substitute structure 
that is selected in such a manner as to give results that 
correspond with reasonable accuracy to the critical 
stresses in the actual structure. 

The art consists in selecting the appropriate idealized 
substitute structure. Primary consideration should 
thereby be given to various factors listed in reference | 
as follows: (1) behavior of thin sheet in compression 
and as a shear web; (2) effects of shear lag; (3) effects 
of doors, windows, and similar cutouts. 


1 The Behavior of Thin Sheet in Compression and as a 

Shear Web 

The behavior of thin sheet in compression has been 
studied extensively elsewhere, and the problem is sum- 
marized in reference | as follows: 

When a compressive load is applied to a stiffened 
panel in the direction of the stiffeners, the thin plate 
will be subjected to a stress of varying intensity. The 
intensity of this stress will tend to be a maximum at 
the stiffeners and a minimum midway between the 
stiffeners. 

For convenience it is assumed that the thin sheet, 
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which is under variable stress, can be replaced by ef- 
fective widths of sheet w cooperating with each stiffener 
and subjected to the same stress as the stiffener, the 
the sheet being considered ineffective—i.e., 
as carrying no load. The equation for computing the 
increments of effective width w can be expressed as 


rest of 


follows: 
w= CtVE/f, 

where ¢ is the thickness of the sheet and f, is the stress 
in the stiffener. Values for the constant C for different 
types of sheet-stiffener combinations using aluminum 
alloys are given in reference 1. 

The behavior of thin sheet acting as a shear web is 
discussed in references 2 and 3. Consideration of this 
effect is a refinement beyond the scope of this paper. 


2 Effects of Shear Lag 

In order to evaluate properly the effects of shear lag, 
it is necessary to have an understanding of the funda- 
mentals of the phenomenon. 

Consider the equilibrium of an element of a cantilever 
prismatic beam, as shown in Fig. 1, with a concentrated 
load P at the free end. The bending moment will be: 


Face A: Px = M 
P (x + Ax) M+AM 


II 


Face B: 


Element AB has been redrawn in Fig. 2 and divided 
into layers of infinitesimal thickness dy. The layers 
have been separated to show their reciprocal effects. 


REFERENCE POINT FOR EQUATION Bf (ZM-0) 
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Normal stress distribution (plane): 
Face A: og = My/I=ce 
Face B: og = (M+ AM)y/I =o + Ao 


Shear stress distribution (constant along x): 
7 = 0.5P(b? — y?)/I 


Equations of equilibrium for Element (1): 


YX = —(o, + Aoijady + osady + 7,a(Ax) = 0 (I) 
SY = +AP, — AP, =0 (Note: TAP=P) (II) 
<M = 7,a(Ax)dy/2 — AP:(Ax) = 0 (IIT) 
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Eqs. (I1) and (III) are of no further interest in the 
matter to be investigated. 

In Eq. (I), the original component of the normal] load 
o,ady is balanced between Face A and Face B longi 
tudinally through Element (1) in pure tension. 

The shear load 7,a( Ax) due to the change in norma! 
load goes first into the neighboring Element (2), from 
there as a component of the shear load 72a( Ax) into the 
next Element (3), and so on through the beam cross 
section, until it is finally balanced by the shear load 
7;a( Ax) in Element (6). The shear load 72,a( Ax) goes 
through Elements (3) and (4) before it is finally bal 
anced in Element (5). 

In general the change in normal load in a given ele 
ment is balanced by an equal change in normal load of 
opposite sign in an element of equal and opposite static 
moment about the neutral axis. Two such elements 
must therefore be connected in such a manner that an 
exchange of shear can take place. 

In the narrow rectangular beam for which the normal 
stress distribution is plane, the shear load path between 
flanges is straight and sufficiently narrow to prevent a 
dissipation or redistribution of stresses over the con 
stant width a. If, on the other hand, the flow of shear 
between flange elements is restricted and the internal 
distribution of longitudinal shear loads along the 
width a to the flange elements is disturbed, a deviation 
from the plane stress distribution is to be expected. 

Consider, for example, the bending stresses in a 
channel as shown in Fig. 3. The shear load in Element 
(1) is ultimately balanced according to the conception 
advanced in the preceding paragraph by the shear load 
in Element (7). However, the shear load path between 
these two flange elements is not straight but leads 
around the corners through Elements (2), (3), (4), (5) 
and (6). Furthermore, all longitudinal shear loads 
between the flanges (i.e., between Elements (1), (2), 
(3) and Elements ), (6), (7)) are balanced only 
through a thin web, represented by Element (4). This 
causes a concentration of shear loads in the thin web 


(i 
(5 


and constitutes a considerable restriction in the flow of 


shear as compared to the narrow rectangle. The 
normal stress distribution will therefore not be 
plane. 
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In the subsequent investigation of the stress distribu- 
tion in a box-beam-type wing, the problem is reduced 
to the study of a channel, since by placing the open ends 
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of two channels together a symmetric box beam is 
formed. 

Consider, therefore, the channel shown in Fig. 4 as 
the left half of a cantilever box beam. The applied load 
is assumed as acting at the shear center, causing pure 
bending. There will therefore, by definition of the shear 
center, be no rotational deflections of the beam, i.e., 
Element d will not deflect laterally and requires no 
support from the right half of the box beam, which is 
therefore justifiably cut away. 

The channel is constructed of thin sheet and stiffened 
by the identical stiffeners a, b, c, d, which are equally 
spaced along the flanges. 

Under the applied load shown in Fig. 4 the upper 
flange is in tension, which is applied to the flange by the 
web along stringer a. In Fig. 5 the longitudinal shear 
load applied by the web to the entire upper flange 
(ADEH) is designated by P,. 
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In Fig. 6 the flange ADEH is shown deflected under 
load, each element being isolated, with the vertical 
loads in equilibrium. Due to the continuous build-up 
of load along the fixed edge EFGH, stringer a shows the 
greatest deflection, stringers }, c and d being next in the 


order named. Hence, from Fig. 6: 


AX’ > Be’ > CC’ > DP! 


Since by Hooke’s Law the stress is proportional to 
the deflection, it follows that 


Ga > oD s To os Od 


where o is the normal stress in the stringer designated 
by the subscript. 
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Since in Fig. 6 


Si = So+ Ri + QO» 


it follows that the load S, applied to Panel I is larger 
than the load S. applied to Panel IT. The shear de- 
flection of Panel I will therefore be greater than the 
shear deflection of Panel II, and the slope of A’B’ with 
respect to AB will be greater than the slope of B’C’ 
with respect to BC. Hence, the difference in deflec- 
tion of stringers a and 0) will be greater than the dif- 
ference in deflection of stringers } and c. 

The same argument applies naturally to Panels II 
and III. Hence, in terms of linear deflections it follows 
that 


(AA’ — BB’) > (BB’ — CC’) > (CC’ — DD’) 


and, as above, in terms of stresses: 


(Oy _ Op) > (op — Ge) > (o, oa) 


This equation for the character of the normal stress 
distribution in the flange is plotted in Fig. 7. 

If area is added to the capstrip stringer a making 
this stringer relatively stiffer, its initial deflection AA’ 
will be smaller. Hence, the initial deflection of Panel 
I will also be smaller, as will that of Panel II and 
so on through the entire flange. Since all deflec- 
tions are smaller, the stress in the entire flange will also 
be lower, as indicated in broken lines in Fig. 7. It 
should be noted, however, that notwithstanding the 
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lower stress in stringer a, the load Q; now carried by 
this stringer due to added area is larger than before, 
since the loads in all other stringers with unchanged 
areas are less, but the total load carried at each longi- 
tudinal station of the flange has not changed. 

If area is added to any other stringer, the load Q 
taken out in that stringer is correspondingly greater 
and will cause a larger change in the angle of the de- 
formed edge A’B’C’D’ at that stringer—the angle 
on the left being larger, and the one on the right being 
smaller, than in the case of identical stringers. 

A change in stringer spacing will affect the stiffness 
of the panel—the closer the stringer spacing, the stiffer 
the panel and the smaller the drop in load carried by 
the neighboring stringer. 

For the complete box beam, which, as stated above, is 
obtained by placing the open ends of two channels to- 
gether, the normal stress distribution for bending about 
a horizontal and vertical axis will be as sketched in 
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By calling the normal stress in the capstrips 100, the 
curves show the ratio of stress in any flange element 
relative ‘to the capstrips and are therefore termed 
“efficiency curves.’ Numerical values taken from these 
curves are known as “efficiency factors.”’ 

In the analysis of a box-beam wing, such efficiency 
curves giving the character of the stress distribution 
can be sketched for various cross sections along the 
wing span. Spanwise plots of the efficiency factors 
along each stringer, which must be fair curves, will 
coordinate the various cross sections. 
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The efficiency factors E obtained from these curves 
are included in the standard bending stress computa- 
tions as follows: 

Moments of inertia: 


I, = SE,Ay’ 

I, = 2E,Ax? 
Normal stresses: 

o, = E,M,y/I, 

o, = E,M,x/I1, 


3 Effects of Doors, Windows, and Similar Cutouts 


By carrying the above reasoning a little further, 
efficiency curves may be drawn for sections through 
cutouts. The efficiency of discontinued stringers is 
naturally zero at the cut, while the efficiency of the 
stringers directly adjacent to a cutout will be consider 
ably higher because of stress concentrations at the 


edges. The spanwise curves along the individual 
stringers must be faired carefully. However, the 
estimation of stress concentrations is beyond the 


scope of this paper. 


THE STRESS DISTRIBUTION IN A BOX BEAM WING AS A 
Case STUDY 


The application of efficiency factors is demonstrated 
for the case of the cantilever wing beam sketched in 
Fig. 9, which has a semispan of 21 ft. and a root chord 
of 6 ft. The wing consists of a center section and two 
outer panels. There are two engine nacelles in the 
center section. The material is Alcoa 24S-T through- 
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Fic. 9. Test specimen. 

The center section is straight, of constant cross sec- 
tion (N.A.C.A. Airfoil 23015), having no taper, sweep- 
back, or dihedral. The structure is semimonocoque, 
having two webs. The upper surface between the webs 
consists of two extruded plates with integral stiffeners. 
The plates are butt-jointed along the center of the box 
by a tee stringer. 

The bending stresses in a critical loading condition 
are calculated for a representative cross section with 














ALLOWANCE FOR SHEAR LAG 


and without use of efficiency factors, and the computa- 
tions are checked by actual stress measurements on a 
full-size test specimen. 


CALCULATION OF STRESSES 


Applied Loads 


The applied loads on the wing at a station 43 in. out- 
board of the airplane centerline for a hydraulic test 
pressure of 700 Ibs. per sq. in. are given in Fig. 10 in 
terms of the following components: 


L = 8,944 lbs. 

D = —2,730 lbs. 

M, = 1,926,600 in.Ibs. 
Mp = —442,600 in.lbs. 


The pitching moment is not of interest and is neglected. 


APPLIED SPANWISE LOAD DISTRIBUTION 
| — \ AND BENDING MOMENT CURVES AT 700 FSi 
HYDRAULIC TEST PRESSURE 
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Fic. 10. 


Selection of Critical Stresses 


The closed type of box beam being studied here has 
a large included area, which is inherently rigid against 
torque; the torsional stresses are small and of no conse- 
quence. 

The shear stresses must usually be carefully checked, 
but in small and medium-size airplanes the shear webs 
are inclined to be stronger than calculated by routine 
conservative engineering methods. that 
may cause trouble here are likely to be limited to the 
design of efficient methods of connecting the shear webs 
to the flanges. 


Problems 
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Since the beam-bending stress is an inverse function 
of the moment of inertia and the height of the wing 
beam is limited, bending strength is dependent on the 
amount of area in the flanges. Sufficient bending 
strength can be provided only at direct cost in weight. 
The chord bending will directly add to the beam-bend- 
ing stress. 

The bending stresses are therefore the primary de- 
sign criterion and comprise the subject of the present 
limited investigation. 


Estimation of Effective Areas 


The trailing edge of the wing aft of the rear web is 
built of relatively thin sheet and will therefore carry 
negligible loads. The leading edge is extensively cut 
out so that major loads at station 43 cannot occur. All 
material forward of the front web and aft of the rear 
web will therefore be considered ineffective, and only 
the flange material between the two webs is assumed 
to carry load. } 

Because of its large thickness, the top cover is as- 
sumed to be totally effective in compression across its 
full width. The lower flange is assumed to be fully 
effective in tension across its full width. 

The effective areas of the webs in compression are 
computed by the formula taken from reference 1 in 
heading No. J. The distribution of effective material 
is shown in Fig. 11. 
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Distribution of effective material in bending. 


Estimation of the Effects of Shear Lag 


Allowance for the effects of shear lag is made in the 
efficiency curves, Fig. 12. These curves were plotted 
after consideration of the arguments presented above, 
and they are claimed to be essentially correct in char- 
acter. The choice of the numerical values is a matter 
of judgment and experience, such as the choice of the 
effective material. 


Computation of Stresses 


Normal bending stresses are computed from the 
values given in Tables 1 and 2 both by the standard 
method and by the suggested method allowing for the 
effects of shear lag. The results of the calculation are 
summarized in Table 3 and plotted in Fig. 13. 
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Fic. 12. Estimated efficiency curves for the effective cross 
section shown in Fig. 11 


TABLE | 
Geometric Values and Efficiency Factors 
Area A, 

Element Sq.In. 2", tn:* 7, %” £y Bm By, % 
(a) 0.759 —16.0 1.2 100 L100 
(b) 0.490 —13.4 8.3 100 97 
(c) 0.400 —11.4 7.95 100 93 
(d) 0.780 — 8.4 8.6 100 88 
(e) 0.400 — 5.4 8.15 100 S4 
(f) 0.760 — 2.4 8.6 100 &2 
(g) 0. 234 0 8.0 100 81 
(h) 0.760 2.4 8.5 100 82 
(j) 0.400 5.4 7.8 100 84 
(k) 0.780 8.4 8.16 100 88 
(m) 0.400 11.4 7.43 100 93 
(n) 0.490 13.4 7.76 LOO 97 
(p) 0.655 15.95 6.70 100 100 
(q) 0.474 16.1 —4.47 100 100 
(r) 1.224 11.5 —5.5 100 94 
(s) 0.152 5.4 —5.0 100 89 
(t) 1.120 0 —5.7 100 86 
(u) 0.152 — § —4.85 100 89 
(v) 1.224 —11.5 —5.1 100 94 
(w) 0. 500 — 16.05 —3.68 100 100 

Total 12.154 


* x’ and y’ are measured from the reference axes, Fig. 11. 
EXPERIMENTAL DETERMINATION OF STRESSES 


Description of the Test Specimen 


The actual bending stresses are determined on a full- 
size test model, consisting of a complete center section 
and a left outer panel. The right outer panel is replaced 
by a jury structure. The stresses are measured on the 
left side of the center section with Huggenberger 
Tensometers located on centerlines of the longitudinal 
members at and between the front and rear web cap- 
strips on both upper and lower surfaces (Fig. 9). 

The model is supported by a jig at the wing-fuselage 
fittings 
Description of the Loading Equipment 

The test loads are applied by hydraulic jacks. The 
applied jack-bending moments, Fig. 10, are corrected 
for the effects of wing dead weight. 
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Fic. 13. Stress distribution over the wing chord at 43 in 
outboard from airplane centerline and at 700 lbs. per sq. in. hy 
draulic test pressure 

TABLE 2 


Moments of Inertia and Applied Moment Vectors 


Standard Suggested 


Item Method Method 
Moments of inertia about r 498.3 $58.2 
reference axes ly, 1,368.4 1,368.4 
Polar moment of inertia k = Axy,* —31.6 —31.4 
Angle of principal axes at —2°4’ —i° oi" 
Moments of inertia about I. 497.4 457.5 
principal axes # 1,368.3 1,370.9 
Applied moment vectors about M,’ 1,910,700 1,912,400 
principal axes Wp’ —511,3800 —504,200 
* Where x1 = x — (ZAx)/(ZA) and y, = y — (TAy)/(ZA) 
} From tan 2a = 2k/(I,, — I;,). 
TABLE 3 
Calculated Stresses 
Normal Normal 
Stress, Stress, 
Standard Suggested 
Fiber “. 50.* +. Ty.” Method Method 
(a) 16.0 —4.4 — 22,900 — 24,280 
(c) 11.4 —5.0 — 23,470 — 23,580 
(e) 5.4 —5.4 — 22,770 — 20,980 
(g) 0 —5.5 — 21,100 — 18,600 
(p) —16.1 —5.1 — 13,580 — 15,400 
(q) —16.5 7.3 34,170 36,550 
(s) 5.9 8.3 34,110 33,000 
(u) 4.9 8.5 30,870 29,800 
(w) 15.5 7.6 23,400 26,120 


* x and y are measured from the principal axes to the centroids 
of elements as designated in Fig. 11. 


The effects of the nacelle dead weight are substituted 
by hand-operated hydraulic ‘“‘nacelle jacks,’ which 
apply down (tension) loads at the center of gravity of 
the nacelles through Chatillon traction dynamometers. 
The nacelle jack loads are applied to the test specimen 
through the wing-nacelle fittings. 











ALLOWANCE 


The “‘lift jacks’ act vertically upward, normal to the 
floor, and are distributed along the wing span so as to 
give the desired vertical bending moment /;. The 
chordwise location determines the applied pitching 
moment. The “drag jacks’’ act upward, normal to 
the floor, on the bell cranks, which pull forward on 
cables attached to the wing. Fig. 14 shows a general 
view of the setup. The equipment is described in 
detail in reference 4. 





Lift and drag jacks operate on a common oil line. 
Since the friction in these specially designed jacks is 
negligible, the oil pressure in all jacks is the same; 
hence, the load applied by the jacks is proportional to 
the jack size and, in the case of the drag jacks, the bell- 
crank ratio. 

This setup insures a simultaneous application of 

The setup is simple to operate and 
The loading is accurate and easy to 


loads by all jacks. 
easy to oversee. 

control. The specimen is readily accessible at all times 
and under all loadings. 


and most of the danger and confusion resulting from 


There are no inertia effects, 
secondary failures is eliminated. 


Interpretation of Test Data 


The test data include stress readings taken at seven 
Huggenberger stations as located in Fig. 9 for a series 
of 46 runs throughout a loading range of 0 to 725 Ibs. 
per sq.in. hydraulic test pressure. From these data, 
straight-line load-stress curves are plotted for each 


Huggenberger station. For convenience, the applied 
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load is given in terms of the hydraulic pressure in the 
jacks. 

The actual stress distribution, curve Fig. 13, is plotted 
for an arbitrarily chosen applied test pressure of 700 Ibs. 
per sq.in. by taking the stresses for each Huggenberger 
station from the corresponding load-stress curve. 


CONCLUSIONS 


The curves of stress distribution calculated accord 
ing to the suggested method follow the shape of the 
measured stress distribution curves considerably better 
than the curves calculated by the standard method. By 
a better choice of the efficiency curves, the exact shape 
of the measured curve could be reproduced. 

The actual values given by the new curve are, in 
general, too high, indicating that the estimation of 
effective areas, Fig. 11, is conservative for this new 
method of calculation. 

The measured values in Fig. 13 for the upper rear 
and lower front corners are higher than the values cal- 
culated the inclination of the 
calculated curves is due to the effect of bending in 
chordwise the actual due to the 
chordwise bending moment Wp are smaller than cal- 
This indicates that the estimation of effective 


in both cases. Since 


direction, stresses 
culated. 
area with respect to chordwise bending is conservative 
and the effect of the trailing edge is not negligible. 

It has therefore been demonstrated that curves of 
approximately identical numerical values and shapes 
as the measured curves can be calculated. 

It is expected that a man of good judgment will be 
able, with some experience, to estimate eifective areas 
and efficiency curves with reasonable accuracy, so that 
this method can be of considerable value in the practical 
design of efficient metal aircraft structures. 

It is, however, highly desirable that reliable methods 
of estimating effective areas and efficiency curves be 


developed. 
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Auditory Sensitivity Under Conditions of 
Anoxia: A Study of Speech Intelligibility’ 


C. P. SEITZ{t anp G. M. SMITH? 
College of the City of New York 


INTRODUCTION 


ee STUDIES CONCERNING the effects of oxy- 
gen deprivation on the organism have been under- 
taken. These studies (especially those of Heymens, 
Cobb, Gildea, McFarland and Evans, Seitz, and Seitz 
and Rosenthal) have implied that the nervous system 
is especially sensitive to lack of oxygen. The response 
of the eye has commonly served as the major indicator 
of decreased neural efficiency. There is reason to be- 
lieve that the ear is also affected by oxygen deprivation 
in a similar manner. A survey of the literature, how- 
ever, indicates that the observed effect of oxygen de- 
privation on the auditory threshold for pure tones is 
almost insignificant except at a relatively high altitude 
20,000 ft., 9.53 per cent oxygen) and this, in spite of 
the fact that Wever and Bray have found marked reduc- 
tion in the cochlear response following a restriction of 
the blood supply by pressure in the arteries. 

In view of these facts and in view of the observed 
variability of the threshold for pure tones at high alti- 
tudes attributable to the distractions induced by an 
adjusting physiology (heart, breathing, etc.), it was 
throught advisable to work with sounds of speech com- 
plexity in the readily audible range. 


PROBLEM 


The purpose of the present study is to inquire into 
the effect of oxygen deprivation on the ability to per- 
ceive speech sounds in the easily audible range of 
stimulus intensities. It is obvious in dealing with this 
problem that two major variables are involved: (1) 
anoxia and (2) pressure. Pressure affects not only the 
availability of oxygen for the organism but also the ear 
and communication system in a mechanical sense. 
This exploratory study is restricted to oxygen depriva- 
tion uncomplicated by pressure changes. 


METHOD AND MATERIALS 


The method of observing the effect of oxygen de- 
privation was to note changes in the subjects’ ability to 
perceive and report words in common speech. In order 
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to obviate the possibility of the experimenter’s attitude 
influencing his production of the speech stimuli, the 
test words were recorded. The recordings were made 
on 16-in. acetate discs with high-fidelity equipment. 
The stimuli were selected from test materials used by 
the Bell Telephone Laboratories* covering most of the 
sounds that occur in common speech. There were seven 
series of approximately 70 monosyllabic words each, 
for both the sea level and the altitude runs. The follow- 
ing combination were employed. 


Test No. 


Description 
1 Single words 
5 Words in groups of 2 
2 Words in groups of 3 
9) 4 Py . 
3 Words in groups of 4 
4 Numbers in groups of 4 


Words in groups of 2 plus range signals 
Words in groups of 2 plus motor noises and 


range 


signals 


The series were distributed so as to minimize the time 
factor, which is important because the effects of oxygen 
deprivation increase with prolonged exposure. 

All of the experiments were carried out in the oxygen 
chamber of the College of the City of New York. The 
chamber is a solid room 8 ft. by 8 ft. by 7 ft. with 
walls approximately 8 in. thick making it relatively 
soundproof. An air-conditioning unit maintained the 
room at comfortable temperature and humidity. An 
altitude of 18,500 ft. (10.2 per cent oxygen) was simu- 
lated by running nitrogen into the chamber and diluting 
the air. This altitude was selected in order to have a 
condition that would be likely to yield unambiguous re- 
sults. Though it is common practice to use oxygen 
above 10,000 or 12,000 ft., a pilot flying at 40,000 to 
45,000 ft. breathing 100 per cent oxygen has a physi- 
ology roughly corresponding to that of a pilot at 16,000 
to 18,000 ft. without oxygen. 

Samples of the chamber air were taken after the as- 
cent and subsequently analyzed on the Haldane-Hen- 
derson-Bailey gas analysis apparatus. None of the 
samples deviated more than the equivalent of 500 ft. 
from the desired altitude. 


SUBJECTS 


Sixteen male college students free from known sys- 
temic defects as determined by medical examination 
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served as subjects. The age range was 17 to 25 with 


the median at 20 years. 


PROCEDURE 


The subjects were examined individually under two 
separate conditions—sea level and altitude. Before the 
experiment proper began, brief instructions were given 
relative to the chamber and its use so as to allay any 
apprehension. The subject was then comfortably 
seated and a pair of Brush Communication Type Head- 
phones* adjusted to the head. An instruction record 
with sample stimuli was then presented and the method 
of response was indicated—simple verbal repetition. 
The experimenter remained in the chamber with the 
subject to reassure him and also to record the responses 
on a check list. During the high-altitude run the ex- 
perimenter was supplied with oxygen through a nasal 
catheter. The ascent took 27 min.; 15 min. were 
spent at altitude to allow for physiologic adjustment 
before any data were collected. Sea level and altitude 
runs were distributed so that each came first half of 
the time. The two forms of the test were also dis- 
tributed so as to balance any inequality in difficulty. 


RESULTS 


Table 1 gives results for each of the 16 subjects on 
tests 1, 2, 3, and 5 (tests without background noises) 

* The response curve of these phones cuts off sharply below 
1,000 e.p.s. and is relatively flat (+10 decibel) from 1,000 to 
Measurements were made at Stevens Institute of 
Technology. The transmission system as a whole was sufficiently 
free from gross distortions to make possible an easy identification 


8,000 c.p.s. 


of the speakers’ voices. 
TABLE | 


The Effect of Altitude on Errors—Tests 1, 2, 3, 5, Combined 








Alt./S.L. 


Subject Sea Level Altitude 
1 29 errors* 32 1.107 
2 22 46 2.09F 
3 44 59 1.34 
4 16 39 2.44 
5 55 41 0.76T 
6 38 21 0.55t 
7 46 67 1.46 
8 52 54 1.04 
9 50 76 1.52 
10 56 77 1.38f 
11 49 43 0.887 
12 35 71 2.03 
13 78 63 0.81t 
14 55 61 1.11 
15 46 87 1.89 
16 52 60 1.15t 
Mean 45.2 56.1 1.35 
Ratio M an./Msg . 1.24 


“ Total number of test items = 280 for each subject (192 
consonants, 88 vowels). 
| Easier form of test taken at altitude; 


€etrors were made at altitude in four of these eight cases. 


in spite of this more 
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in terms of total errors. The total number of items 
was 280, of which 192 tested consonants and 88 vowels. 
It will be seen that in 12 out of 16 cases more errors 
were made at altitude than at sea level. The mean of 
the error ratios is 1.35, indicating a mean increase in 
errors of 35 per cent. The ratio of the means, a more 
stable value, f shows an increase in errors of 24 per cent. 
The individual differences in the error ratios are con- 
siderable, ranging from 2.44 to 0.55. Only four of these 
ratios are less than unity. It is interesting to observe 
that in all of these cases the easier form of the tests 
was taken at high altitude. Undoubtedly, some of the 
variability is a function of this unequal difficulty in the 
two forms of the tests. This inequality does not distort 
the results, however, since half of the subjects took the 
easier form at sea level and half at altitude. It 


reasonable to suppose that a considerable part of the 


is 


variability is due to individual differences in tolerance 
for low oxygen. 

The data of Table 1, expressed in errors, become more 
meaningful if expressed in terms of articulation values 
for vowels, consonants, and standard syllables as in 
Table 2. In doing this the procedure of the Bell Labora- 
tories was followed, where articulation value is defined 
as the per cent of stimuli correctly perceived. The con- 
sonant and vowel articulation were determined directly 
from our error data. The syllable articulations were 
calculated from the Bell Laboratory formula’ derived 
= | 


empirically from extensive observations, S 
(1 — VC?)?-9, 

+ This ratio, obtained by dividing the sum of the errors at alti 
tude by the sum of the errors at sea level, is less affected by ex- 
treme individual variations that enter into the calculation of the 
mean of the error ratios. 


TABLE 2 


Values for Consonants, Vowels, and Standard 


Syllables—Tests 1, 2, 3, 5, Combined 


Articulation 


— —Sea Level— - — Altitude—————. 

Sub- Con- Con- 

ject sonant Vowel Syllable* sonant Vowel Syllable* 
1 85.0% 100.0% 67.0% 85.0% 96.5% 65.0% 
2 91.0 94.5 75.0 77.0 97.5 54.0 
3 77.5 99.0 6.0 70.5 97.5 45.0 
4 92.5 97.5 + 80.0 79.5 100.0 59.0 
5 71.5 100.0 47.5 79.0 99.0 8.0 
6 82.0 96.5 61.0 90.0 97.5 75.5 
7 76.0 100.0 54.0 69.0 92.0 40.5 
8 73.5 99.0 50.0 73.0 97.5 48.5 
9 75.0 98.0 51.0 66.0 87.5 35.0 
10 71.0 100.0 46.5 64.0 91.0 34.0 
11 76.0 98.0 53.0 78.0 100.0 57.0 
12 82.0 100.0 63.0 69.0 86.0 38.0 
13 62.0 94.0 33.0 68.0 98.0 42.0 
14 73.0 95.0 47.5 70.0 79.0 35.5 
15 77.0 99.0 55.0 57.0 94.0 28.0 
16 74.0 98.0 50.5 69.0 100.0 44.0 

Mean 77.5 98.0 55.5 72.8 94.6 47.5 
P.E. x7 1.25 0.3 1.98 1.13 0.49 3.33 
*S=1-—-—(1 — VC). 
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As might be expected, consonants are less readily per- 
ceived than vowels at both sea level and altitude. The 
mean syllable articulation at sea level is 55.5 per cent; 
at high altitude it falls to 47.5 per cent. The reliability 
of the difference between these means was computed by 
Fisher's method for small correlated samples. The / 
value was 2.7, with a P value of 0.0082, indicating that 
the probability is less than 1 in 100 that the observed 
difference was due to chance. 

Though the effect of altitude on syllable articulation 
may not seem impressive, an examination of Fig. 1 re- 
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SYLLABLE ARTICULATION 
Fic. 1. Di crete word intelligibility vs. syllable articulation 
(Reproduced from Bell Telephone Laboratory Bulletin B-436, 
with the permission of the author 


produced from the report of Fletcher and Steinberg? 
will show that a loss in intelligibility of this magnitude 
This figure 


may be critical for word intelligibility. 
shows the relationship between syllable articulation 
and single word intelligibility. A critical point is 
reached in the neighborhood of 30 per cent syllable 
articulation, where a slight change in syllable articula- 
tion corresponds to a much larger change in word in- 


telligibility. No communication system gives perfect 
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intelligibility under conditions of flight because of 
noise, vibration, and low pressure. Any additional fac 
tors such as anoxia or poor speech production may re 
duce articulation to the critical point. 

The results based upon our two tests with background 
noises are mentioned with some diffidence because it 
was felt that the experimental procedure that lies back 
of them was based on inadequate equipment. However 
the procedure did have exploratory value and the results 
suggest further lines for investigation. There was a 
marked drop in consonant, vowel, and syllable arti 
culation as might have been anticipated, the mean 
syllable articulation values being reduced to about 20 
per cent. The altitude effect was small but consistent 
with the results for the tests without background noises 

It is noteworthy that no errors were made in the 
This is 


to be expected since numbers can be identified from 


number series at either sea level or altitude. 


the vowel sounds alone and it will be recalled that vowel 
articulation was 98 per cent at sea level and 95 per cent 
at altitude. This is of considerable importance in work 

ing out codes such as are now being adopted in military 
and civil aviation in order to comply with restrictions 


imposed by wartime conditions. 
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